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APPLICATIONS OF MODERN HYDRODYNAMICS TO AERONAUTICS. 

By L. PRANDTL. 


PART I. 

FUNDAMENTAL CONCEPTS AND THE MOST IMPORTANT THEOREMS. 

1. All actual fluids show internal friction (viscosity), yet the fofces due to viscosity, with 
the dimensions and velocities ordinarily occurring in practice, are so very small in comparison 
with the forces due to inertia, for water as well as for air, that we seem justified, as a first ap- 
proximation, in entirely neglecting viscosity. Since the consideration of viscosity in the 
mathematical treatment of the problem introduces difficulties which have so far been overcome 
only in a few specially simple cases, we are forced to neglect entirely internal friction unless we 
wish to do without the mathematical treatment. 

We must now ask how far this is allowable for actual fluids, and how far not. A closer 
examination shows us that for the interior of the fluid we can immediately apply our knowl- 
edge of the motion of a nonviscous fluid, but that care must be taken in considering the layers 
of the fluid in the immediate neighborhood of solid bodies. Friction between fluid and solid 
body never comes into consideration in the fields of application to be treated here, because it 
is established by reliable experiments that fluids like water and air never slide on the surface 
of the body; what happens is, the final fluid layer immediately in contact with the body is 
attached to it (is at rest relative to it), and all the friction of fluids with solid bodies is therefore 
an internal friction of the fluid. Theory and experiment agree in indicating that the transition 
from the velocity of the body to that of the stream in such a case takes place in a thin layer of 
the fluid, which is so much the thinner, the less the viscosity. In this layer, which we call the 
boundary layer, the forces due to viscosity are of the same order of magnitude as the forces due 
to inertia, as may be seen without difficulty. 1 It is therefore important to prove that, however 
small the viscosity is, there are always in a boundary layer on the surface of the body forces 
due to viscosity (reckoned per unit volume) which are of the same order of magnitude as those 
due to inertia. Closer investigation concerning this shows that under certain conditions there 
may occur a reversal of flow in the boundary layer, and as a consequence a stopping of the fluid 
in the layer which is set in rotation by the viscous forces, so that, further on, the whole flow is 
changed owing to the formation of vortices. The analysis of the phenomena which lead to the 
formation of vortices shows that it takes place where the fluid experiences a retardation of flow 
along the body. The retardation in softie cases must reach a certain finite amount so that a 
reverse flow arises. Such retardation of flow occurs regularly in the rear of blunt bodies; there- 
fore vortices are formed there very soon after the flow begins, and consequently the results 
which are furnished by the theory of nonviscous flow can not be applied. On the other hand, 
in the rear of very tapering bodies the retardations are often so small that there is no noticeable 
formation of vortices. The principal successful results of hydrodynamics apply to this case. 
Since it is these tapering bodies which offer specially small resistance and which, therefore, 
have found special consideration in aeronautics under similar applications, the theory can be 
made useful exactly for those bodies which are of most technical interest. 

1 From this consideration one can calculate the approximate thickness of the boundary layer for each special case. 
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For the considerations which follow we obtain from what has gone before the result that 
in the interior of the fluid its viscosity, if it is small, has no essential influence, but that for 
layers of the fluid in immediate contact with solid bodies exceptions to the laws of a nonviscous 
fluid must be allowable. We shall try to formulate these exceptions so as to be, as far as possi- 
ble, in agreement with the facts of experiment. 

2. A further remark must be made concerning the effect of the compressibility of the 
fluid upon the character of the flow in the case of the motion of solid bodies in the fluid. All 
actual fluids are compressible. In order to compress a volume of air by 1 per cent, a pressure 
of about one one-hundredth of an atmosphere is needed. In the case of water, to produce an 
equal change in volume, a pressure of 200 atmospheres is required; the difference therefore is 
very great. With water it is nearly always allowable to neglect the changes in volume arising 
from the pressure differences due to the motions, and therefore to treat it as absolutely incom- 
pressible. But also in the case of motions in air we can ignore the compressibility so long as 
the pressure differences caused by the motion are sufficiently small. Consideration of compressi- 
bility in the mathematical treatment of flow phenomena introduces such great difficulties that 
we will quietly neglect volume changes of several per cent, and in the calculations air will be 
looked upon as incompressible. A compression of 3 per cent, for instance, occurs in front of a 
body which is being moved with a velocity of about 80 m./sec. It is seen, then, that it appears 
allowable to neglect the compressibility in the ordinary applications to technical aeronautics. 
Only with the blades of the air screw do essentially greater velocities occur, and in this case the 
influence of the compressibility is to be expected and has already been observed. The motion 
of a body with great velocity has been investigated up to the present, only along general lines. 
It appears that if the velocity of motion exceeds that of sound for the fluid, the phenomena are 
changed entirely, but that up close to this velocity the flow is approximately of the same char- 
acter as in an incompressible fluid. 

3. We shall concern ourselves in what follows only with a nonviscous and incompressible 
fluid, about which we have learned that it will furnish an approximation sufficient for our 
applications, with the reservations made. Such a fluid is also called “ the ideal fluid.” 

What are the properties of such an ideal fluid ? I do not consider it here my task to develop 
and to prove all of them, since the theorems of classical hydrodynamics are contained in all 
textbooks on the subject and may be studied there. I propose to state in what follows, for 
the benefit of those readers who have not yet studied hydrodynamics, the most important 
principles and theorems which will be needed for further developments, in such a manner that 
these developments may be grasped. 1 ask these readers, therefore, simply to believe the 
theorems which I shall state until they have the time to study the subject in some textbook 
on hydrodynamics. 

The principal method of description of problems in hydrodynamics consists in expressing in 
formulas as functions of space and time the velocity of flow, given by its three rectangular com- 
ponents, u, v, u and in addition the fluid pressure p. The condition of flow is evidently com- 
pletely known if u, v , w t and p are given as functions of x } y, 2 , and t, since then u , v, w , and p 
can be calculated for any arbitrarily selected point and for every instant of time. The direc- 
tion of flow is defined by the ratios of u, v, and w; the magnitude of the velocity is \ u 2 -f v 2 -f u> 2 . 
The “streamlines” will be obtained if lines are drawn which coincide with the direction of 
flow at all points where they touch, which can be accomplished mathematically by an inte- 
gration. If the flow described by the formulas is to be that caused by a definite body, then 
at those points in space, which at any instant form the surface of the body, the components of 
the fluid velocity normal to this surface must coincide with the corresponding components 
of the velocity of the body. In this way the condition is expressed that neither does the fluid 
penetrate into the body nor is there any gap between it and the fluid. If the body is at rest 
in a stream, the normal components of the velocity at its surface must be zero; that is, the flow 
must be tangential to the surface, which in this case therefore is formed of stream lines. 
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4. In a stationary flow— that is, in a flow which docs not change with the time, in which 
then every new fluid particle, when it replaces another particle in front of it, assumes its velocity , 
both in magnitude and in direction and also the same pressure there is, for the fluid particles 
lying on the same stream line, a very remarkable relation between the magnitude of the velocity, 
designated here by V, and the pressure, the so-called Bcrnouilli ccpiation 

p t V 1 — const. (1) 

( P is the density of the fluid, i. e., the mass of a unit volume). This relation is at once appli- 
cable to the case of a body moving uniformly and in a straight line in a fluid at rest, for we are 
always at liberty to use for our discussions any reference system having a uniform motion in a 
straight line. If we make the Velocity of the reference system coincide with that of the body, 
then the body is at rest with reference to it, and the flow around it is stationary. If now V 
is the velocity of the bodv relative to the stationary air, the bitter will have in the new refer- 
ence system the velocity V upon the body (a man on an airplane in flight makes observations 
in terms of such a reference system, and feels the motion of flight as “wind ). 

The flow of incident air is divided at a blunt body, as shown in figure 1. At the point A 
the flow comes completely to rest, and then is again set in motion in opposite directions, tan- 
gential to the surface of the body. We learn from equation 
(1) that at such a point, which we shall call a “rest-point,” 

the pressure must be greater by 2 than in the undisturbed 

fluid. We shall call the magnitude of this pressure, of which 
we shall make frequent use ; the “dynamical pressure,” and 
shall designate it by q. An open end of a tube facing the 
stream produces a rest point of a similar kind, and there arises 
in the interior of the tube, as very careful experiments have 
shown, the exact dynamical pressure, so that this principle 
can ho used for the measurement of the velocity, and is in 
fact much used. The dynamical pressure is also well suited 
to express the laws of air resistance. It is known that this resistance is proportional to 

the square of the velocity and to the density of the medium; but q -^V 2 ) so the law of air 
resistance may also be expressed by the formula 

W^c.F.q ( 2 ) 

where F is the area of the surface and c is a pure number. With this mode of expression it 
appears very clearly that the force called the “drag” is equal to surface times piessuie differ- 
ence (the formula has the same form as the one for the piston force in a steam engine). This 
mode of stating the relation has .been introduced in Germany and Austria and has proved use- 
ful. The air-resistance coefficients then become twice as large as the “absolute” coefficients 
previously used. 

Since V 2 can not become less than zero, an increase of pressure greater than q can not, by 
equation (1), occur. For diminution of pressure, however, no definite limit can be set. In 
the case of flow past convex surfaces marked increases of velocity of flow occur and in connection 
with them diminutions of pressure which frequently amount to 3 q and more. 

5. A series of typical properties of motion of nonviscous fluids may be deduced in a useful 
manner from the following theorem, which is due to Lord Kelvin. Before the theorem itself 
is stated, two concepts must be defined. 1. The circulation: Consider the line integral of the 
velocity / V cos (F, ds). ds t which is formed exactly like the line integral of a force, which is 
called “the work of the force.” The amount of this line integral, taken over a path which 
returns on itself is called the circulation of the flow. 2. The fluid line: By this is meant a line 
which is always formed of the same fluid particles, which therefore shares in the motion of the 
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fluid. The theorem of Lord Kelvin is: In a nonviscous fluid the circulation along every fluid 
line remains unchanged as time goes on. But the following must he added: 

(1) The case may arise that a fluid line is intersected by a solid body moving in the fluid. 
If this occurs, the theorem ceases to apply. As an example I mention the case in which one 
pushes a flat plate into a fluid at rest, and then by means of the plate exerts a pressure on the 
fluid. By this a circulation arises which will remain if afterwards the plate is quickly withdrawn 
in its own plane. See figure 2. 

(2) In order that the theorem may apply, we must exclude mass forces of such a character 
that work is furnished by them along a path which returns on itself. Such forces do not ordi- 
narily arise and need not be taken into account here, where we are concerned regularly only 
with gravity. 

(3) The fluid must be homogeneous, i. c., of the same density at all points. We can easily 
see that in the case of non uniform density circulation can arise of itself in the course of time 
if we think of the natural ascent of heated air in the midst of cold air. The circulation increases 
continuously along a line which passes upward in the warm air and returns downward in the 
cold air. 



Frequently the case arises that the fluid at the beginning is at rest or in absolutely uniform 
motion, so that the circulation for every imaginable closed line in the fluid is zero. Our theorem 
then says that for every closed line that can arise from one of the originally closed lines the 
circulation remains zero, in which we must make exception, as mentioned above, of those lines 
which are cut by bodies. If the line integral along every closed line is zero, the line integral 
for an open curve from a definite point 0 to an arbitrary point P is independent of the selection 
of the line along which the integral is taken (if this were not so, and if the 
integrals along two lines from 0 to P were different, it is evident that the 
line integral along the closed curve OPO would not be zero, which contra- 
dicts our premise). The line integral along the line OP depends, therefore, 
since we will consider once for all the point 0 as a fixed one, only on the coordi- 
nates of the point P, or, expressed differently, it is a function of these coor- 
dinates. From analogy with corresponding considerations in the case of 
fields of force, this line integral is called the “ velocity potential,” and the 
particular kind of motion in which such a potential exists is called a “ potential motion.” As 
follows immediately from the meaning of line integrals, the component of the velocity in a 
definite direction is the derivative of the potential in this direction. If the line-element is 
perpendicular to the resultant velocity, the increase of the potential equals zero, i. e., the sur- 
faces of constant potential are everywhere normal to the velocity of flow. The velocity itself 
is called the gradient of the potential. The velocity components u, v f w are connected with the 
potential 4> by the following equations: 


Fig. 2 — Production of c ir- 
culation by introduc- 
lion arid withdrawal of 
Hat plate. 


d<f> d<{> d(f> 

u = ~ 1 ■> 

ox oy Oz 


(3) 


The fact that the flow takes place without any change in volume is expressed by stating that 
as much flows out of every element of volume as flows in. This leads to the equation 


<)u , Ov t Ow _ 

< .. i — () 

Ox d y dz 


(4) 


In the case of potential flow we therefore have 


d 2 4> d 2 0 d 2 0 
dx 2 Oy 2 dz 2 


(4a) 


as the condition for flow without change in volume. All functions (x, y, z, f ), which satisfy 
this last equation, represent possible forms of flow. This representation of a flow is specially 
convenient for calculations, since by it the entire flow is given by means of the one function 4>. 
The most valuable property of the representations is, though, that the sum of two, or of as 
many as one desires, functions 4 J , each of which satisfies equation (4a), also satisfies this equation, 
and therefore represents a possible type of flow (“ superposition of flows”). 
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6. Another concept can be derived from the circulation, which is convenient for many con- 
siderations, viz, that of rotation. The component of the rotation with reference to any axis 
is obtained if the circulation is taken around an elementary surface of unit area in a plane 
perpendicular to the axis. Expressed more exactly, such a rotation component is the ratio of 
the circulation around the edge of any such infinitesimal surface to the area of the surface. The 
total rotation is a vector and is obtained from the rotation components for three mutually per- 
pendicular axes. In the case that the fluid rotates like a rigid body, the rotation thus defined 
comes out as twice the angular velocity of the rigid body. If we take a rectangular system of 
axes and consider the rotations with reference to the separate axes, we find that the rotation can 
also be expressed as the geometrical sum of the angular velocities with reference to the three axes. 

The statement that in the case of a potential motion the circulation is zero for every 
closed fluid line can now be expressed by saying the rotation in it is always zero. The theorem 
that the circulation, if it is zero, remains zero under the conditions mentioned, can also now 
be expressed by saying that, if these conditions are satisfied in a fluid in which there is no 
rotation, rotation can never arise. An irrotational fluid motion, therefore, always remains 
irrotational. In this, however, the following exceptions arc to be noted: If the fluid is divided 
owing to bodies being present in it, the theorem under consideration does not apply to the 
fluid layer in which the divided flow r reunites, not only in the case of figure 2 but also in the 


case of stationary phenomena as in figure 3 , 
since in this case a closed fluid line drawn in 
front of the body can not be transformed into 
a fluid line that intersects the region where the 
fluid streams come together. Figure 3 shows 
four successive shapes of such a fluid line. This 
region is, besides, filled with fluid particles which 
have come very close to the body. We are 
therefore led to the conclusion from the stand- 
point of a fluid with very small but not entirely 
vanishing viscosity that the appearance of vor- 
tices at the points of reunion of the flow in the 
rear of the body does not contradict the laws of 
rotation £, 77, f arc expressed as follows by mcai 



body. 

hydrodynamics. The three components of the 
is of the velocity components u, v , w. 


dw dv _du dw, dv _du 
n ~dz dx 9 * dy 

If the velocity components arc derived from a potential, as shown in equation (2), the rotation 

. . . , . n . d 2 * _ d 2 < f> 

components, according to equation ( 5 ) vanish identically, since 

7 . Very remarkable theorems hold for the rotation, which were discovered by v. Helmholtz 
and stated in his famous work on vortex motions. Concerning the geometrical properties of the 
rotation the following must be said: 

At all points of the lluid where rotation exists the direction of the resultant rotation axes 
can be indicated, and lines can also be drawn whose directions coincide everywhere with these 
axes, just as the stream lines are drawn so as to coincide with the directions of the velocity. 
These lines will be called, following Helmholtz, “vortex lines/ The vortex lines through the 
points of a small closed curve form a tube called a “vortex tube.” It is an immediate con- 
sequence of the geometrical idea of rotation as deduced above that through the entire extent 
of a vortex tube its strength - i. e., the circulation around the boundary of the tube is constant. 
It is seen, in fact, that on geometrical grounds the space distribution of rotation quite in e- 
pendently of the special properties of the velocity field from which it is deduced is of the same 
nature as the space distribution of the velocities in an incompressible fluid. Consequently a 
vortex tube, just like a stream line in an incompressible fluid, can not end anywhere m the 
interior of the fluid; and the strength of the vortex, exactly like the quantity of fluid passing 
per second through the tube of stream lines, has at one and the same instant the same va ue 
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throughout the vortex tube. If Lord Kelvin’s theorem is now applied to the closed fluid line 
which forms the edge of a small element of the surface of a vortex tube, the circulation along it 
is zero, since the surface inclosed is parallel to the rotation axis at that point. Since the circula- 
tion can not change with the time, it follows that the element of surface at all later times will 
also be part of the surface of a vortex tube. If we picture the entire bounding surface of a vortex 
tube as made up of such elementary surfaces, it is evident that, since as the motion continues 
this relation remains unchanged, the particles of the fluid which at any one time have formed 
the boundary of a vortex tube will continue to form its boundary. From the consideration 
of the circulation along a closed line inclosing the vortex tube, we see that this circulation — i. e., 
the strength of our vortex tube — has the same value at all times. Thus we have obtained the 
theorems of Helmholtz, which now can be expressed as follows, calling the contents of a vortex 
tube a “vortex filament “The particles of a fluid which at any instant belong to a vortex 
filament always remain in it; the strength of a vortex filament throughout its extent and for 
all time has the same value. ” From this follows, among other things, that if a portion of the 
filament is stretched, say, to double its length, and thereby its cross section made one-half as 
great, then the rotation is doubled, because the strength of the vortex, the product of the rota- 
tion and the cross section, must remain the same. We arrive, therefore, at the result that the 
vector expressing the rotation is changed in magnitude and direction exactly as the distance 
between two neighboring particles on the axis of the filament is changed. 

8. From the way the strengths of vortices have been defined it follows for a space filled 
with any arbitrary vortex filaments, as a consequence of a known theorem of Stokes, that 
the circulation around any closed line is equal to the algebraic sum of the vortex strengths 
of all the filaments which cross a surface having the closed line as its boundary. If this closed 
line is in any way continuously changed so that filaments are thereby cut, then evidently the 
circulation is changed according to the extent of the strengths of the vortices which are cut. 
Conversely we may conclude from the circumstance that the circulation around a closed line 
(which naturally can not be a fluid line) is changed by a definite amount by a certain displace- 
ment, that by the displacement vortex strength of this amount will be cut, or expressed differ- 
ently, that the surface passed over by the closed line in its displacement is traversed by vortex 
filaments whose strengths add up algebraically to the amount of the change in the circulation. 

The theorems concerning vortex motion are specially important because in many cases 
it is easier to make a statement as to the shape of the vortex filaments than as to the shape of 
the stream lines, and because there is a mode of calculation by means of which the velocity 
at any point of the space may be determined from a knowledge of the distribution of the rota- 
tion. This formula, so important for us, must now r be discussed. If F is the strength of a 
thin vortex filament and ds an element of its medial line, and if, further, r is the distance from 
the vortex element to a point P at which the velocity is to be calculated, finally if a is the angle 
between ds and r, then the amount of the velocity due to the vortex element is 


dv 


F ds sin a, 
4 7 r r 3 ; 


( 6 ) 


the direction of this contribution to the velocity is perpendicular to the plane of ds and r. The 
total velocity at the point P is obtained if the contributions of all the vortex elements present 
in the space are added. The law 7 for this calculation agrees then exactly with that of Biot- 
Savart, by the help of which the magnetic field due to an electric current is calculated. Vor- 
tex filaments correspond in it to the electric currents, and the vector of the velocity to the 
vector of the magnetic field. 

As an example w T e may take an infinitely long straight vortex filament. The contributions 
to the velocity at a point P are all in the same direction, and the total velocity can be deter- 
mined by a simple integration of equation (G). Therefore this velocity is 

+ 'X 

r C ds - sin a 

V = 4xJ 
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As seen by figure 4, s = h ctg a, and by differentiation, ds= — ^~da. Further r = 


so that 


T 

, = r r 

4 7 r hj 


sura 


sin ada= — 4 ^-^[cos a]* = 2 ^ 


(6 a) 


This result could be deduced in a simpler manner from the concept of circulation if we were 
to use the theorem, already proved, that the circulation for any closed line coincides with 
the vortex strength of the filaments which are inclosed by it. The circulation for every closed 
line which goes once around a single filament must therefore coincide with its strength. If 
the velocity at a point of a circle of radius h around our straight filament equals v then this 

F 

circulation equals “path times velocity ” = 2 irh v, whence immediately follows v = (j^j l m The 

more exact investigation of this velocity field shows that for every point outside the filament 
(and the formula applies only to such points) the rotation is zero, so that in fact we are treat- 
ing the case of a velocity distribution in which only along the axis does rotation prevail, at 
all other points rotation is not present. 

For a finite portion of a straight vortex filament the preceding calculation gives the value 


v = 4tt 7i,( cos ai ~ ^ os 


(6b) 


This formula may be applied only for a scries of portions of vortices which together give an 
infinite or a closed line. The velocity field of a single portion of a filament would require 
rotation also outside the filament, in the sense that from the end of 
the portion of the filament vortex lines spread out in all the space 
and then all return together at the beginning of the portion. In the 
case of a line that has no ends this external rotation is removed, 
since one end always coincides with the beginning of another portion 
of equal strength, and rotation is present only where it is predicated 
in the calculation. 


h 


-4- 


J 

^ds - 


FIG. F— Yelocif y-lield due to infinite 
rectilinear vortex. 


9. If one wishes to represent the flow around solid bodies in a fluid, one can in many cases 
proceed by imagining the place of the solid bodies taken by the fluid, in the interior of which 
disturbances of flow (singularities) are introduced, by which the flow is so altered that the 
boundaries of the bodies become streamline surfaces. For such hypothetical constructions 
in the interior of the space actually occupied by the body, one can assume, for instance, any 
suitably selected vortices, which, however, since they are only imaginary, need not obey the 
laws of Helmholtz. As we shall sec later, such imaginary vortices can be the seat of lifting 
forces. Sources and sinks also, i. e., points where fluid continuously appears, or disappears, 
offer a useful method for constructions of this kind. While vortex filaments can actually 
occur in the fluid, such sources and sinks may be assumed only in that part of the space which 
actually is occupied by the body, since they represent a phenomenon which can not be realized. 
A contradiction of the law of the conservation of matter is avoided, however, if there are assumed 
to be inside the body both sources and sinks, of equal strengths, so that the fluid produced by 
the sources is taken back again by the sinks. 

The method of sources and sinks will be described in greater detail when certain practical 
problems are discussed; but at this point, to make the matter clearer, the distribution of veloci- 
ties in the case of a source may be described. It is very simple, the flow takes place out from 
the source uniformly on all sides in the direction of the radii. Let us describe around the point 
source a concentric spherical surface, then, if the fluid output per second is Q f the velocity at 
the surface is 

iirr 2 ’ 


(7) 
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the velocity therefore decreases inversely proportional to the square of the distance, 
flow is a potential one, the potential comes out (as line-integral along the radius) 

Q 

= const. — . 

4tt r 


The 

(7a) 


If a uniform velocity toward the right of the whole fluid mass is superimposed on this 
velocity distribution — while the point source remains stationary — then a flow is obtained 
which, at a considerable distance from the source, is in straight lines from left to right. The 
fluid coming out of the source is therefore pressed toward the right (see fig. 5); it fills, at some 
distance from the source, a cylinder whose diameter may be determined easily. If V is the 
velocity of the uniform flow, the radius r of the cylinder is given by the condition Q^irr 2 - V. 
All that is necessary now is to assume on the axis of the source further to the right a sink of 
the same strength as the source for the whole mass of fluid from the source to vanish in this, 
and the flow closes up behind the sink again exactly as it opened out in front of the source. 
In this way we obtain the flow around an elongated body with blunt ends. 


10. The special case when in a fluid flow the phenomena in all planes which are parallel 
to a given plane coincide absolutely plays an important role both practically and theoretically. 

If the lines which connect the corresponding 
points of the different planes are perpendicular to 
the planes, and all the streamlines are plane 
curves which lie entirely in one of those planeS, 
we speak of a uniplanar flow. The flow around 
a strut whose axis is perpendicular to the direc- 
tion of the wind is an example of such a motion. 

The mathematical treatment of plane poten- 
tial flow of the ideal fluid has been worked out 
specially completely more than any other prob- 
lem in hydrodynamics. This is due to the fact 
that with the help of the complex quantities 
(x + iy, where i = \ — 1 , is called .the imaginary 
unit) there can be deduced from every analytic 
function a case of flow of this type which is incom- 
pressible and irrotational. Every real function, 4> (x, y) and 'k (x, y), which satisfies the relation 

4>-fi4'==/(x + iy), (8) 

where f is any analytic function, is the potential of such a flow. This can be seen from these 
considerations: Let x + iy be put=z, where z is now a lt complex number.” Differentiate equa- 
tion (8) first with reference to x and then with reference to y, thus giving 

.d* _df dz = df 
dx ^ % dx dz dx dz 



Fig. ■'). — Superposition of uniform flow and that caused by a source. 


d4> . c N' d/dz _ ,df_ .04* _ d * 

by 1 by dz by ~ % dz " 7 dx by 

In these the real parts on the two sides of the equations must he equal and the imaginary 
parts also. If 4> is selected as the potential, the velocity components u and v are given by 


_d4>_ M 
U dx by ’ V by dx 


(9) 


If now we write the expressions 
4* and then of they become 


bu bv 
dx by 


(continuity) and 


bv 

dx 


(rotation) first in terms of 
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d U ^ dv _ d 2 4> <44 <44 (44 ' 

dx dy dx 2 dy 2 dydx dxdy 
= o 

dv_du_d 2 <t> d 2 4 _ <44 <44 

dx dy dydx dxdy~ dx 2 dy 2 

=--o 


(10) 


It is seen therefore that not only is the motion irrotational (as is self-evident since there is 

a potential), hut it is also continuous. The relation ^ ^ = 0 besides corresponds exactly 

to our equation (4a). Since it is satisfied also by 4, this can also be used as potential. 

The function 4', however, has, with reference to the flow deduced by using 4 as potential, a 
special individual meaning. From equation (8) we can easily deduce that the lines 4 — const, 
are parallel to the velocity; therefore, in other words, they are streamlines. In fact if we put 


d>V = 'J dx 4- f ^dy = o } t hen = 
dy J 5 dx 


d4 . 
dx° 


d 4 


dx 

d4 

dy 


v 

u 


which expresses the fact, of parallelism. The lines 4 — const, are therefore perpendicular to the 
lines 4 = const. If we draw families of lines, 4= const. and 4= const. for values of 4 and 4 
which differ from each other by the same small amount, it follows from the easily derived 


equation (74 4 7(74 


dx + idy ) that the two bundles form a square network; from wliie 


h fol- 


lows that the diagonal curves of the network again form an orthogonal and in fact a square 
network. This fact can be used practically in drawing such families of curves, because an error 
jn the drawing can be recognized by the eye in the wrong shape of the network of diagonal 
curves and so can he improved. With a little practice fairly good accuracy may be obtained 
by simply using the eye. Naturally there are also mathematical methods for further improve- 
ment of such networks of curves. The function 4, which is called the “stream function,” 
has another special meaning. If we consider two streamlines 4 - 4, and 4 = 4 2 , the quantity 
of fluid which flows between the two streamlines in a unit of time in a region of uniplanar flow 
of thickness 1 equals 4 2 — 4,. In fact if we consider the flow through a plane perpendicular to 
the .Y-axis, this quantity is 

r*\h r*y* pyi 

Q= I udy - / d *jdy I d'l' 'K - 'f', . 

U v\ U yi J j/i 


The numerical value of the stream function coincides therefore with the quantity of fluid which 
flows between the point j, y and the streamline 4 = o. 

As an example let the function p 

4 4 i4 = A (x+ iy ) n 


be discussed briefly. It is simplest in general to ask first about the streamline 4 = o. As is 
well known, if a transformation is made from rectangular coordinates to polar ones r, <p, (x*F iyY 
= r n (cos ny> 4 i sin n<p). The imaginary part of this expression is ir n sin ntp. This is to be 
put equal to i4. 4 = o therefore gives sin n<p = o, i. e., n<p~o, 7 r, 2 tt, etc. The streamlines 4 = o 


are therefore straight lines through the origin of coordinates, which make an angle a = ^ with 

each other, the flow is therefore the potential flow between two plane walls making the angle 
a with each other. The other streamlines satisfy the equation r n sin ny> = const. The veloci- 
ties can be obtained by differentiation, e. g., with reference to x: 


0<P .<54 

dx 4 1 dx “ U ~ 1 V = "1 71 ( X 4 b ? /) 


n— 1 


= A nr° 1 { co.s (a-1) 4 i sin (/i-l )<^ } 
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For r = o this expression becomes zero or infinite, according as n is greater or less than 11, i. e., 
according as the angle ct is less or greater than ir(= 180°). Figures 6 and 7 give the streamlines 

7T 3 9 

for a = ^ = 45° and 9 7 t = 270°, corresponding to n — 4 and Jn the case of figure 7 the velocity, 

as just explained, becomes infinite at the corner. It would be expected that in the case of 
the actual flow some effect due to friction would enter. In fact there are observed at such 
corners, at the beginning of the motion, great velocities, and immediately thereafter the for- 
mation of vortices, by which the motion is so changed that the velocity 
at the corner becomes finite. 

It must also be noted that with on cquatiorf 

p+iq = <p{x + iy) (11) 

the x-y plane can be mapped upon the p-q plane, since to every pair 
of values x;y a pair of values p,q corresponds, to every point of the x-y 
plane corresponds a point of the p-q plane, and therefore also to every 
element of a line or to every curve in the former plane a linear element 
and a curve in the latter plane. The transformation keeps all angles 
unchanged, i. e., corresponding lines intersect in both figures at the same 
angle. 

By inverting the function <p of equation (11) we can write 

x + iy = x (p+iq) 

and therefore deduce from equation (8) that 

<T> + 'i'k =f [x (p 4- iq) ] = F (p 4- iq) (12) 

and v are connected therefore with p and q by an equation of the type of equation (8), and 
hence, in the p-q plane, are potential and stream functions of a flow, and further of that flow 
which arises from the transformation of the 4\ ^ network in the x-y planp into the p-q plane. 

This is a powerful method used to obtain by transformation from a known simple flow 
new types of flow for other given boundaries. Applications of this will be given in section 11 



Fig. (i. — riiiplutjar Mow he- 
t wet’ll plane walls making 
an angle a— 4;'° with each 
other. 


11. The discussion of the principles of the hydrodynamics of nonviscous fluids to be 
applied by us may be stopped here. I add but one considera- 
tion, which has reference to a very useful theorem for obtaining 
the forces in fluid motion, namely the so-called “momentum theo- 
rem for stationary motions.” 

We have to apply to fluid motion the theorem of general 
mechanics, which states that the rate of change with the time 
of the linear momentum is equal to the resultant of all the ex- 
ternal forces. To do this, consider a definite portion of the 
fluid separated from the rest of the fluid by a closed surface. 

This surface may, in accordance with the spirit of the theorem, 
he considered as a “fluid surface,” i. e., made up always of 
the same fluid particles. We must now state in a formula the 
change of the momentum of the fluid within the surface. If, as 
we shall assume, the flow is stationary, then after a time (it every fluid particle in the interior 
will be replaced by another, which has the same velocity as had the former. On the boundary, 
however, owing to its displacement, mass will pass out at the side where the fluid is approaching, 
and a corresponding mass will enter on the side away from which the flow takes place. If dS 
is the area of an element of surface, and v n the component of the velocity in the direction of 
the outward drawn normal at this element, then at this point dm^pdS . v n dt. If we wish 
to derive the component of the “ impulse” — defined as the time rate of the change of momen- 
tum— for any direction s, the contribution to it of the element of surface is 



Fig. 7. — Uniplanar flow around plane 
walls making an angle 270° wit h each 
other. 


dJ t = v^-pdS 


v„ v. 


(13) 
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With this formula wc have made the transition from the fluid surface to a corresponding solid 
“ control surface.” 

The external forces are compounded of the fluid pressures on the control surface and the 
forces which are exercised on the fluid by any solid bodies which may be inside of the control 
surface. If we call the latter P, we obtain the equation 

2 *WJ V • cos (r?, s) . dS +p ( 14 ) 

for the s component of the momentum theorem. The surface integrals are to be taken over the 
entire closed control surface. The impulse integral can be limited to the exit side, if for every 
velocity on that side the velocity v is known with which the same particle arrives at the 
approach side. Then in equation (13) dJ is to be replaced by 

dJ - dJ r = ( v 8 — v s ') ^ = pdSv n (v B — % 0 (13a) 

The applications given in Part II will furnish illustrations of the theorem. 
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PART II. 

APPLICATIONS. 


A. DISTRIBUTION OF PRESSURE ON AIRSHIP BODIES. 

12. Thu first application of hydrodynamical theory to be tested by experiment in the 
Gottingen Laboratory referred to the distribution of pressure over the surface of models of 
airships. We can construct mathematically the flow for any number of varieties of sectional 
forms of botlies of revohition of this kind if we place along an axis parallel with the direction 
of the air current any suitable distribution of sources and sinks, taking care that the total 
strength of the sources and sinks are the same. According to the intensity of the uniform 
motion which is superimposed upon the flow from the sources, we obtain from the same system 
of sources and sinks bodies of different thicknesses. In order to obtain the smoothest possible 

shapes, the sources and sinks are generally distributed contin- 
uously along the axis, although single-point sources are allowable. 

In the case of continuously distributed'sourccs and sinks the 
method of procedure is briefly this: The abscissas of the single 
sources are denoted by £, the intensity of the source per unit of 
length by /(£), in which positive values of /(£) denote sources, nega- 
tive values sinks. The condition that makes the stream from the 
sources self-contained is expressed by the equation 





Fit;. *.'■ -Explanation of qnanUtu^ 
in rak-uhilion of streamline shapes 
Abscissae: Position of source. 
Ordinates: Intensity of source per 
length, i. c. f($). 


(‘ /({)(/£ = 0 . 
wo 


By simply adding the potentials due to the single elementary 
sources /(£) i. e., in this case by integrating them, the total flow due to the sources will be 
given by the potential defined by the following formula 


(x,y) = 


if 




05 ) 


in which v (/--£) 2 + and y is the perpendicular distance from the axis of the point for 
which the potential is calculated, x is the abscissa along the axis measured from the same 
origin as £. (See fig. 8.) There must be added to this potential that due to the uniform flow 
with the velocity V, viz, <f>,= Vx. The total potential is then ( l> = <1^ f <L; and therefore the veloc- 


r)<t) rl(t> 

V+ n 1 and the sidewise (radial) velocity is v= .. = 1 

J dy dy 


ity parallel to the axis is ^ 

In order to calculate the streamlines one could perform an integration of the direction 
given by u and v. These lines are obtained more conveniently, in this case also, by means of 
the stream function. (See sec. 10.) In the case of flow symmetrical with reference to the axis, 
such as is here discussed, one can take as stream function the quantity of fluid flowing inside 
the circle drawn through the point x, y in a plane perpendicular to the axis and having its center 
172 
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on the axis. The amount of fluid delivered by the sources which lie up the stream is purposely 
deducted from this. It is not difficult to see that all points of the X—Y plane, through whose 
parallel circles the same amount of fluid flows per second — after deduction of the sources — 
must lie on one and the same streamline, for evidently there is no flow, cither in or out, through 
the surface formed by the streamlines drawn through the points of any one parallel circle 
(since the flow is along the surface) ; therefore the quantity of fluid flowing within this surface is 
constant, so far as it is not increased by the sources. From the meaning of the stream function, 
to determine which the velocity must be integrated over a surface, it follows that the stream 
function of a flow due to two or more causes is at every point the sum of the stream functions 
of the several partial flows. For a continuous distribution of sources therefore the stream 
function \p is obtained by an integration exactly as was the potential. According to our 
premise the surface of the body is designated simply by the value \p = o. The formulas arc 
obtained as follows: 

The flow from a simple source through a circle passing through a point lying to the right of 
the source is, writing r = ^x 2 -\- y 2 , 






From this, in accordance with what has been said, the quantity Q must be subtracted, so that 




<2 = _ 2 ( 1+ r) 

For points lying to the left of the source we obtain from the integral 

Qt 


(10) 


"T 


?K) 


(17) 


which coincides with formula (1G); this holds, then, everywhere. 

For the assumed continuous distribution of sources we obtain 

*■— iJX'+V)* 

in which £) 2 + ?/ 2 - To this stream function of the sources must now be added that due 

to the paralhl flow ^ (18) 

Putting the total stream function 4q 4 4U — equal to zero, gives the equation of the surface 
of the body around which the flow takes place. Putting gives any other streamline. 

It is evident that, with the same distribution of sources, a wholes group 
of body surfaces can be obtained, depending upon the choice of the ratio 
of the intensity of the sources to the strength of the parallel flow. 

The determination is best made practically by graphical methods, for 
Instance, by laying off the curves x = const, in a system of coordinates 
consisting of y and which can be obtained at once from a calculation 
by tables for the stream function Sfq. If we intersect these curves by 
parabolas corresponding to the equation - -ty — W?/ 2 — C, we obtain at 
once a contour (for C=o), or some external or internal streamline (for 
C>o or C<o). The parabola maybe drawn upon transparent paper, and 
then by displacing the parabola along the 'F axis we can at once obtain 
from figure 9 the values of y corresponding to any x. 

In this manner a former colleague of mine, who unfortunately fell immediately at the 
beginning of the war, Dr. G. Fuhrmann, calculated the shapes of bodies corresponding to a 
series of source distributions, and on the one hand he determined the distribution of pressure 
over the surface of these bodies by means of the Bernouilli equation (see sec. 4) 



Fig. 9. — The curves are fur 
different values of x = 
const. 


P = Po+o( V 2 - (u 3 + D 3 )}‘ 


(19) 


i The velocities u and v may be obtained from the potential, but also from the stream function 


i/ ; f6r u - 


i (i\i> , l d’p 

„ J and r«= - 

2i ry r iy 2wt/ 


2i ry dx 
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and on the other he constructed models according to these drawings and measured the pressure 
distribution over them when placed in a wind tunnel. The agreement was altogether surpris- 
ingly good, and this success gave us the stimulus to seek further relations between theoretical 
hydronamics and practical aeronautics. The work of Fuhrmann was published in Jahr-b. 
der Motorluftschiff-Studien Gesellsch., Volume V, 1911-12 (Springer, Berlin), and contains 
a large number of illustrations. Four of the models investigated are shown here. The upper 
halves of figures 10 to 13 show the streamlines for a reference system at rest with reference to 
the undisturbed air, the lower halves the streamlines for a reference system attached to the 
body. The distribution of the source intensities is indicated on the axis. The pressure dis- 
tributions are shown in figures 14 to 17. The calculated pressure distributions are indicated 
by the lines which are drawn full, the individual observed pressures by tiny circles. 2 

It is seen that the agreement is very complete; at the rear end, however, there appears a 
characteristic deviation in all cases, since the theoretical pressure distribution reaches the full 
dynamical pressure at the point where the flow reunites again, while actually this rise in pres- 
sure, owing to the influence of the layer of air retarded by friction, remains close to the surface. 

As is well known there is no resistance for the theoretical flow in a nonviscous fluid. The 
actual drag consists of two parts, one resulting from all the normal forces (pressures) acting 
on the surface of the body, the other from all the tangential forces (friction). The pressure 
resistance, which in this case can be obtained by integration of the pressure distribution over 
the surface of the body, arises in the main from the deviation mentioned at the rear end, and 
is, as is known, very small. Fuhrmanrds calculations gave for these resistances a coefficient, 
with reference to the volume of the body, as shown in the following table: 3 


Model 

1 

II 

III 

IV 




ki 


0.0170 

0.0123 

0.0131 

0. 0145 


This coefficient is obtained from the following formula: 

Drag w l =Jc 1 IJ 2 > 3 q 

where V designates the volume and q the dynamical pressure. 

The total resistance (drag) was obtained for the four models by means of the balance; 
the difference between the two quantities then furnishes the frictional resistance. The total 
drag coefficients were : 3 


Model 


k 


I 

0.0340 


II 

0. 0220 


III | IV 
0. 0240 0. 0248 


With greater values of VL than were then available for us, the resistance coefficients 
become nearly 30 per cent smaller. For purposes of comparison with other cases it may be 
mentioned that the “maximum section' 5 was about 2/5 of U 2 > 3 . The surface was about seven 
times U 23 ; from which can be deduced that the total resistance of the good models was not 
greater than the friction of a plane surface having the same area. The theoretical theorem 
that in the ideal fluid the resistance is zero receives in this a brilliant confirmation by experiment. 

B. THEORY OF LIFT. 

13. The phenomena which give rise to the lift of an aerofoil may be studied in the simplest 
manner in the case of uniplanar motion. (See sec. 10.) Such a uniplanar flow would be ex- 
pected obviously in the case that the wing was unlimited at the sides, therefore was “ infinitely 

1 In the wind tunnel there was a small pressure drop in the direction of Us length. In order to eliminate the effect of this, the pressures toward 
the fore had to be diminished somewhat and those aft somewhat increased. 

* After deduction of the horizontal buoyancy. 
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long,'* and throughout exhibited the same profile and the same angle of attack. In this case all 
the sections will be alike in all respects and each one can be considered as a plane of symmetry. 
The infinitely long wing plays an important part therefore in the considerations of the theoretical 
student. It is not possible to realize it in free air. and marked deviations from the infinitely 
long wing are shown even with very long wings, e. g., those having an aspect ratio of 1 : 10. In 
laboratories, however, the infinitely long wing, or uniplan ar flow, may be secured with good 
approximation, if a wing having a constant profile is placed between plane walls in a wind 
tunnel, the walls running the full height of the air stream. In this case the wing must extend 
close to the walls; there must be no gap through which a sensible amount of air can flow. We will 

now discuss such experiments, and first we shall state the funda- 
mental theory of uniplanar flow. 

Since, as explained in section 4, in a previously undisturbed 
fluid flow, the sum of the static and dynamic pressures is con- 
stant: p- j- ^ F 2 = const., in order to produce lift, for which the 

pressure below the surface must be increased and that above 
diminished, such arrangements must he made as will diminish 
the velocity below the wing and increase it above. The other 
method of producing such pressure differences, namely, by 
fig. is.— Deduction of the Kutia formula, causing a vortex region above the surface placed lvi*e a kite 
unipinnar now around iTinmii* wm K . oblique to the wind, by which a suction is produced, does not 

come under discussion in practical aeronautics owing to the great resistance it sets up. Lanchester 
has already called attention to the fact that this lifting current around the wing arises if there is 
superimposed upon a simple potential flow a circulating flow which on the pressure side runs 
against the main current and on the suction side with it. Kutta (1902) and Joukowski (1906) 
proved, independently of each other, the theorem that the lift for the length l of the wing is 

A = pV VI ( 20 ) 

in which r is the circulation of the superimposed flow. It may be concluded from this formula 
that in a steady fluid flow lift is not possible unless there is motion giving rise to a circulation. 
In uniplanar How in an ideal fluid this lift does not entail 
an} 7 drag. 

The proof of the Kutta- Joukowski formula is generally 
deduced by applying the momentum theorem to a circular 
cylinder of large radius whose axis is the medial line of the 
wing. The circulatory motion, which could he obtained 
numerically close to the wing only by elaborate mathematical 
processes, is reduced at a great distance from the wing to a 
motion which agrees exactly with the flpw around a rectilinear 
vortex filament (see sec. 8), in which, therefore, the single i'i<= so. rmphnmr iuiir»rm flow anumi dr- 
particles describe concentric circles. The velocity around a niI:,r 

circle of radius R is, tlten, v — For an element of surface l. RdB (see fig. 18) the normal 

component of the velocity is V cos 9, the mass flowing through per second dm — plRV cos 8d6. 
If we wish to apply the momentum theorem for the vertical components, i. e., those perpen- 
dicular to the direction of F, then this component of the velocity through the element of surface 
must be taken. This, obviously, is v cos 9 , taken positive if directed downward; the total 
impulse, then, is 

J — jv cos 9dm = plR FvJ** r cos 2 Odd. 

The integral equals ir, and therefore introducing the value of v 

J-l'VVL 
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Since the resulting impulse is directed downward (the upward velocity in front of the wing is 
changed into a downward one behind the wing), this means that the reaction of the fluid against 
the wing is a lift of the wing upward. The amount of the impulse furnishes, as is seen by re- 
ferring to formula (20), only half the lift. The other half comes from the pressure differences 
on the control surfaces. Since, for a sufficiently large II, v can always lie considered small com- 
pared with Vj neglecting £ v 2 , the pressure p is given, according to the crnouilli equation, by 

p ----- p 0 + 2 V* - 2 1 ( v s ' m 02 + v 2 cos 2 ~ (• V sin 0. 

A component of this, obtained by multiplying by sin 9, acts vertically on the surface element 
IRdd. The resulting force I) is, then, 

I) ~ plR Vvj q sin 2 Odd. 

This integral also equals tt, so that here also 

D = J pVTl 

its direction is vertically up. The total lift, then, is 

A =J + I) = pW . 

14. For the more accurate analysis of the flow around wings the complex functions (see 
sec. 10) have been applied with great success, following the procedure of Kutta. Very different 



Fig. 20. — Unipluimr flow around circular 
cylinder considered as a columnar vortex 
of strength r. 



methods have been used. Here we shall calculate only one specially simple case, in which the 
flow will be deduced first around a circular cylinder and then calculated for a wing profile by a 
transformation of the circular cylinder and its flow, using complex functions. 

The flow around a circular cylinder has long been known. If the coordinates in the plane 
of the circle are p and q, and if we write p + iq = t, the potential and stream functions for the 
ordinary symmetrical flow around the circular cylinder are given by the very simple formula 

+ 5 ( 21 ) 

It is easily seen by passing to polar coordinates that, for r = a, \p { =^o, and that therefore the 
circle of radius a is a streamline. Further, for the p axis, ^= 0 , i. e., this is also a streamline. 
The whole flow is that shown in figure 19. To this flow must be added the circulation flow 
expressed by the formula 

$1 + i% = 2 W lo S 1 ' (22) 

r 

which, as shown in figure 20, is simply a flow in concentric circles with the velocity 9 . The 

combination of the two flows, i. e., the flow for the sum of the expressions in equations (21) 
and (22), is shown in figure 21. It is seen that the rest point is moved down an amount de. 
By a suitable choice of the circulation this can be brought to any desired point. 

a i / a *\ / a*\ 

» t+ t = ( r+ f ) cos t?+i ( r- f i sin 8 i log t=- tf+t log r. 

20167—23 12 
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We must now discuss the transformation of this How to a wing profile. For this purpose 
manifold means are possible. The simplest is furnished by a transformation according to the 
equation 

z = x + iy = t + Y' 

By this the circle of diameter AB=^2b in the t plane (as we shall for brevity’s sake call the p f q 
plane) is transformed into a straight line A f B f of the length 4b along the X axis, and concen- 



tric circles around the former become ellipses, the radii become hyperbolas. All the ellipses 
and hyperbolas have their foci at the ends of the straight line, this forming a confocal system. 
Figures 22 and 23 illustrate the transformation. It may be mentioned, in addition, that the 
interior of the circle in figure 22 corresponds to a continuation of the meshwork in figure 23 


through the slit A f B\ whose form agrees with the meshwork as drawn. Any circle through 



Via. 2L— IUust rations of Joukowski sections. 


the points A B is thereby transformed into an arc of 
a circle passed over twice, having an angle subtended 
at the center equal to 4/3. 

Many different results may now be obtained by 
means of this mapping, according to the position 
which the circle, around which the flow takes place 
according to equations (21) and (22), bears to the di- 
ameter AB of the circle of figure 22. If the diameter 
AB is made to coincide with any oblique diameter of 
the circular section of the cylinder, we obtain a flow 
around an oblique plate whose angle of at tack coin- 
cides with the inclination of the line AB. If the di- 


ameter AB is selected somewhat smaller, so that both points lie inside the circle symmetrically 
on the diameter, the flow around ellipses is obtained. If, however, the diameter AB coincides 


with a chord of the circle around which the original flow was, which, for example, may lie below 


the center, the flow around a curved plate forming an arc of a circle is obtained. By selection 
of various points in the interior of the original circle forms of diverse shapes are obtained. The 
recognition of the fact that among these forms very beautiful wunglike profiles may be found we 
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owe to Joukowski. These are obtained if the point B is selected on the boundary of the original 
circle and the point A inside, and somewhat below the diameter through the point B. Figure 
24 gives illustrations of such Joukowski profiles. 

In order that the flow may be like the actual one, in the cases mentioned the circulation 

must always be so chosen that the rear rest 

point coincides with the point B } or, re- - 

spectively, with the point on the original - — ^ - 

circle which lies nearest this point. In ___ 

this case there will be, after mapping on — - — *7 — - — 

the z plane, a smooth flow away from the — " — - — 

trailing edge, as is observed in practice. It 

is therefore seen that the circulation must 

be taken greater according as the angle of ^ — 

attack is greater, which agrees with the ob- 

servation that the lift increases with increas- — 

ing angle of attack. 

The transformation of the flows shown T , , 4 . Ifl 

_ . Eicj. 'j.i. - I ruiNfurumtiun of smiplo potential flow, figure 10. 

in figures 19 to 21 into wing profiles gives 

illustrations of streamlines as shown in figures 25 to 27 — figure 25, simple potential flow; 
figure 2(3, circulation flow; figure 27, the actual flow around a wing obtained by superposition 
of the two previous flows. 

We are, accordingly, by the help of such constructions, in the position of being able to 

S calculate the velocity at every point in the 
neighborhood of the wing profile, and with it 
v the pressure. In particular, the distribution 
\ of pressure over the wing itself may be cal- 
\ culated. 

\ My assistant, Dr. A. Betz, in the year 

I 1914 worked out the pressure distribution for 
1 a Joukowski wing profile, for a series of angles 
/ of attack, and then in a wind tunnel meas- 
J ured the pressure distribution on a hollow 
/ model of such a wing made of sheet metal, 
side walls of the height of the tunnel being 

Fkj. --Trims formation of circulatory flow, figure 20. introduced SO AS to Secure Uniplaiiar flow 

The results of the measurements agreed in a very satisfactory manner with the calculations, 
on ly — a s could be well explained as due to friction— the actual circulation was always slightly 
less than that calculated for the same angle of attack. If the pressure distributions would 


— Tnmsfortnjlion of circulatory flow, figure 20. 


be compared, not for the same angles of 
attack, but for the same amouht of circu- 
lation, the agreement would he noticeably 
better. The pressure distributions 'are 
shown in figures 28 to 30, in which again 
the full curves correspond to the measure- 
ments and the dashes to the calculated 
pressures. Lift and drag for the. wing 
were also obtained by the wind-tunnel bal- 
ance. In order to do this, the middle part 
of the wing was isolated from the side parts, 
which were fastened to the walls of the tun- 



ncl by Carefully designed labyrinths, so that Flo. 27.— Transformation of superposition of the two flows, figure 21. 

within a small range it could move without 

friction. The result of the experiment is shown in figure 31. The theoretical diag is zeio, that 
obtained by measurement is very small for that region where the wing is “good,” but sensibly 
larger for too large and too small angles of attack. The lift is correspondingly in agreement 
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with the theoretical value in the good region, only everywhere somewhat less. The deviations 
of drag as well as of lift are to be explained by the influence of the viscosity of the fluid. The 
agreement on the whole is as good as can be expected from a theory which neglects completely 
the viscosity. 

For the connection between the angle of incidence a and the circulation which results from 
the condition discussed above calculations give the following result for the lift: 

(1) The Kutta theory gives for the thin plane 

I plate the formula 

| A =bt.Trp]i 2 sin a (23) 

T The lift coefficient 6 T a is defined by the equation 


A 1 

f a = j? where q= 9 p V 2 




and therefore 


= 2tt sin a 


(2) For the circularly curved plates having an 
angle of arc 40 subtended at the center (see figure 23) 
we have, according to Kutta, if a is the angle of 
attack of the chord, 


= 2ir 


sin (a 4-0) 


which, for small curvatures, becomes 2?r sin (a-f0); 
j this can be expressed by saying that the lift of the 

circularly curved plate is the same as that of a plane 
nr;. which touches the former at a point three-fourths of 

the distance around the arc from its leading edge. 
U"l\ oc ^ 6 For the Joukowski profiles and for others the 

-j* formulas are less simple, v. Mises showed in 1917 

' ° that the increase of C a with the angle of attack, i. e., 

a * * s greater for all other profiles than for the flat 

j/// plate, and is the greater the thicker the profile. But 

// the differences arc not marked for the profiles occur- 

ring in practice. 

The movement of the center of pressure has 
j / also been investigated theoretically. With the 

|V / plane plate, in the region of small angles, it always 

| / lies at one-fourth of the width of the plate; with 

| j circularly curved thin plates its position for small 

1 1 angles is given by the following law: 

/ tan a 

11G ‘ ,i0 ' X ° 4 tan a 4- tan 0 

Pressure distribution over a Joukowski wing, different angles of . . 

attack. Full lines give results of wind-tunnel tests; dashed m which t IS the chord of tile plate, and X 0 is the dis- 

lmes, calculated values. tance measured from the center of the plate. The 

fact that the movement of the center of pressure in the case of “good” angles of attack of the 
profiles agrees with theory is proved by the agreement of the actual pressure distribution with 
that calculated. In the case of thin plates a less satisfactory agreement as respects pressure 
distribution is to be expected because with them in practice there is a formation of vortices at 
the sharp leading edges, while theory must assume a smooth flow at this edge. 

15. That a circulatory motion is essential for the production of lift of an aerofoil is defi- 
nitely established. The question then is how to reconcile this fact with the proposition that 
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the circulation around a fluid line in a nonviscous fluid remains constant. If, before the motion 
begins, we draw a closed line around the wing, then, so long as everything is at rest, the circula- 
tion certainly is zero. Even when the motion begins, it can not change for this line. The ex- 
planation of why, in spite of this, the wing gains circulation is this: At the first moment of the 
motion there is still no circulation present, the motion takes place approximately according to 
figure 25, there is a flow at high velocity around the 
trailing edge. (See sec. 10.) This motion can not, how- 
ever, continue; there is instantly formed at the trailing 
edge a vortex of increasing intensity, which, in accord- 
ance with the Helmholtz theorem that the vortex is 
always made up of the same fluid particles, remains with 
the fluid as it passes on. (See fig. 32.) The circulation 
around the wing and vortex, taken together, remain:', 
equal to zero; there remains then around the wing a cir- 
culation equal and opposite to that of the vortex which 
has gone off with the current. Therefore vortices will 
be given off until the circulation around the wing is of 
such a strength as to make the fluid flow off smoothly 
from the trailing edge. If by some alteration of the 
angle of attack the condition for smooth flow is ‘dis- 
turbed, vortices are again given off until the circulation 
reaches its new value. These phenomena are com- 
pleted in a comparatively short distance, so the full lift 
is developed very quickly. 

In the pictures of flow around a wing, e. g., figure 27, 
one sees that the air in front of the wing flows upward 
against the reaction of the lift. The consideration of 
momentum has shown that half of the impulse is due to 
the oncoming ascending current. This fact needs some 
further explanation. The best answer is that given by 
Lanchcster, 7 who shows that for the production of lift 
the air mass at any time below the wing must be given 
an acceleration downward. The question he asks is: What kind of a motion arises if for a short 
time the air below the wing is accelerated downward, then the wing is moved forward a bit 
without pressure, then the air is again accelerated, and so on? The space distribution of the 
accelerations is known for the case of a plane plate, infinitely extended at the sides, accelerated 

= ^ from rest; the pattern of the accelcra- 

tion direction is given in figure 33. It 
/ \t \ is seen that above and below the plate 

\ the acceleration is downward, in front 
i of and behind the plate it is upward 
j opposite to the acceleration of the 
/ plate, since the air is escaping from 
/ the plate. Lanchcster asks now about 

^ the velocities which arise from the 

ITr,. n. -ITiidiuM ion of circMilutinn around a v.iriR due to vortices having (roiling original Uniform Velocity relative to 

edRPl the plate owing to the fact that the 

plate, while it gives rise to the accelerations as shown in figure 33, gradually comes nearer the 
air particle considered, passes by it, and finally again moves forward away from it. The pic- 
ture of the velocities and streamlines which Lanchcster obtained in this way and reproduced 
in his book was, independently of him, calculated exactly by Kutta. It is reproduced in figure 
34. It is seen that as the result of the upward accelerations of* the flow away from the wing 

1 Aerodynamics T, § 110-J1K. 
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there is an upward velocity in front* of the plate, a uniform downward acceleration at the plate 
itself due to which the upward velocity is changed into a downward one, and finally behind the 
plate a gradual decrease of the downward velocity on account of the acceleration upward. 

C. THE FINITE WING. 



Fig. ID. — Acceleration 
diagram around an 
infinitely long flat 
plate accelerated at 
right angle to its sur- 
face. 


16. It has been known for a long time that the aspect ratio of an aerofoil had a great effect 
on its properties. One could therefore have expected that, on account of the vanishing of 
pressure at the side edges, the intensity of the lift must decrease toward the edge, so that its 
average value for the same angle of attack must be smaller for small values 
of the aspect ratio than for large ones. But the observed influence of aspect 
ratios is sensibly greater than could be explained in this way. We must 
therefore investigate whether an explanation of this phenomenon can be found, 
if we apply to the finite aerofoil in some proper manner the results which are 
known to hold for uni planar flow'. 

It is easily seen that vortices in the free fluid must here be taken into 
account. For it is certain that circulation is present around the middle of 
the wing, because no lift is possible without circulation. If a closed line 
drawn around the middle of the wing, around which, therefore, there is circulation, is displaced 
sideways over the end of a wing, it will certainly no longer show 7 circulation here when it is 
beyond the w 7 ing. From the theorem that the circulation along a closed line only changes 
if it cuts vortex filaments, and that the amount of the change of the circulation equals the 
sum of the strengths of the vortex filaments cut (see sec, 8), we must conclude that from each 
half of a w r ing vortex filaments whose strengths add up to r must proceed, which are concen- 
trated mainly near the ends of the whig. According to t he Helm- 

holtz theorem we know r further that every vortex produced in the 
fluid continues to move w r ith the same fluid particles. We may 
look upon the velocities produced by the w 7 ing as small compared 
w r ith the flight velocity V, so that as an approximation w 7 e may 
assume that the vortices move aw r ay from the w r ing backwards 
wuth the rectilinear velocity V. (If it is wished, we can also im- 
prove the considerations based upon such an assumption if the 
motion of the vortices of themselves relative to the air is taken into 
account. This will, however, be seen to be unnecessary for practical applications of t lie theory.) 

In order now r to obtain the simplest possible scheme, we shall assume that the lift is uni- 
formly distributed over the w r ing; then the total circulation will arise only at the ends, and 
continue rearwards as free vortices. The velocity field of an infinitely long w 7 ing, as w r e saw, 

w 7 as the same at great distances as that of a rectilinear vortex 
filament instead of the w T ing. We shall assume that the corre- 
sponding statement holds for the finite wnng. We thus obtain, 
for the velocity field around a finite wing, a picture which is 
somewhat crude, it is true, if w r e take for it the velocity distri- 
bution due to a vortex filament of corresponding shape. 

It may be mentioned here that, on account of there being 
the same law 7 s for the velocity field of a vortex filament and the 
magnetic field of an electric current (see sec. 8) 1 , the velocity 
near a finite w T ing can also be investigated numerically by cal- 
culating the direction and intensity of the magnetic field pro- 
duced near an electrical conductor shaped as showm in figure 35 due to an electrical current 
flowing in it. 

The principles for the calculation of this velocity field have been stated in section 8; the 
total velocity is made up out of three partial velocities which are caused by the three rectilinear 
vortex portions. As is seen without difficulty, for the region between the vortices the flow is 
downward, outside it is upw r ard. 



■^ireumlino^ around tm infinitely 
Ion" curved plate. 




CD 


1 0 - 


Fig. 


V finite wing, considered as due to 
vortices replacing the wing. 
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17. This approximation theorem is specially convenient if the conditions at great distances 
from the wing are treated. With its help we can explain how the weight of an airplane is 
transferred to the ground. In order to make the flow satisfy the condition that at the ground 
components of velocity normal to it are impossible, we apply a concept taken from other 
branches of physics and superimpose the condition of an image of the airplane, taking the earth 
as the mirror. On account of symmetry, then, all velocity components normal to the earth s 
surface will vanish. If we use as our system of coordinates one attached to the airplane, we 
have then the case of stationary motion. If we take the X axis in the direction of the span 
of the wing, the Y axis horizontal in the direction of flight and the Z axis vertically down, 
and if u, v y w are the components of the additional velocity due to the vortices, then calling 
p Q the undisturbed pressure and p ' the pres- 
sure difference from p OJ and neglecting the 
weight of the air, Bernouilli’s equation gives us 

p„ + v' + 2 D 2 + {v - F) 2 + to*] = p 0 4 g V 2 

If this is expanded and if n 2 , v 2 } and w 2 are 
neglected as being small of a higher order, there 
remains the simple equation 

p' = P Vv (27) 

For the determination of the pressure dis- 
tribution on the ground we must now calculate 
the value of v. Let us assume the vortices 
run off the wing in an exactly horizontal direc- 
tion (actually, their path inclines downward 
slightly), in which case they do not contribute 
to v. There remains then only the “ transverse 
vortex” of the length l (effective span) and 
the circulation F. We will assume that the 



span of the wing is small in comparison with 
the distance h of the airplane from the ground. 
In that case we can treat the transverse vortex 
as if it were a single vortex element. We ob- 
tain, then — see figure 36 — at a point A, with 



r 


in. 3^.— Application of method of images to airplane flying near the 
ground. 


the coordinates x and y, a velocity perpendicular to the plane ABF , of the amount 


v 1 = r l 


sin a 
4 x R 2 


The image of the airplaite furnishes an equal amount perpendicular to the plane ABF ' 
If is the angle between the plane ABF and the XY plane, then the actual velocity at the 
ground, as far as it is due to the trans'verse vortex, will be the resultant of v t and v 2 . It is there- 

R f h 

fore v = 2v i sin 0, or, if we write sin a = ^ sin i 9 = ^/ ( see fy?* 36) 


rlh 

2 7 r tf 3 


If we take into account the fact that, according to the Kutta-Joukowski formula (20), 
p r Vl = A, equations (27) and (28) iead to the relation 


Ad 1 1 X i' 

If this is integrated over the whole infinite ground surface, it is seen that the resultant 
force due to the pressures on the ground has exactly the amount A. It is thus proved that the 
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pressure distribution due to the circulation motion transfers to the ground exactly the weigh 
of the airplane. The distribution of the pressure, which according to formula (29) is axially 
symmetrical with reference to the foot of the vertical line drawn from the airplane, is shown in 

A 

figure 37. The pressure maximum is p x = ^ Its amount, even for low heights of flight, is 

very small, since the surface over which the pressure is distributed is very large. 

18. Applications of an entirely different kind may be made of the velocity field which belongs 
to the vortex of figure 35. For instance, an estimate may be made as to the magnitude of the 
downward velocity component at any point of the tail surfaces, and in. this manner the influence 
of the wings upon the tail surfaces may be calculated. If in accordance with the Kutta- 
Joukowski formula the lift is written A = pT VI, in which, taking account of the fact that a 



portion of the vortices flow off within the ends of the wing, l , can be taken somewhat less than 
the actual span b , then at a distance d, behind the wing, the velocity component downward is 


in which a = 




(30) 


If the flight velocity is F, this gives for the inclination of the downward sloping air-current 

V) 

tan <p=* y* We proved this relation in the year 1911 and found an approximate agreement 
with observation. 

The principle made use of above has been applied with profit to the calculation of the 
influence of one wing of a biplane upon the other wing and has given a method for the calcula- 
tion of the properties of a biplane from the properties of a single wing as found by experiments. 
The fundamental idea, which is always applied in such calculations, is that, owing to the vortex 
system of one wing, the velocity field near the wing is disturbed, and it is assumed that a wing 
experiences the same lift as in an undisturbed air stream if it cuts the streamlines of the flow 
disturbed by the other wing in the same manner as a monoplane wing cuts the straight stream- 
lines of the undisturbed flow. As is easily seen, the wing profile must in general be slightly 
turned and its curvature slightly altered, as is shown in figures 38 and 39. By the rotation of 
the wing the direction of the resultant air force acting on it is turned through an equal angle. 
If the magnitude of the velocity as well as its direction is also changed, this must be expressed 
by a corresponding change in the resultant air force. 
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As an illustration we will treat briefly the case of a biplane without stagger. The most 
important component of the disturbance velocity w is again the vertical one; in the plane of 
the mean lift lines of the biplane it is affected only by the pair of vortices running off the wings, 
since the transverse vortex of one wing causes only an increase (or decrease) of the velocity of 
flow at the other wing. We are concerned here only with the calculation of that downward 
disturbance velocity due to the vortices from the wing not under investigation, since the other 
vortex system is present with the monoplane and its influence has already been taken into 
account in the experiments on a monoplane. 

The total velocity due to a portion of a vortex proceeding to infinity in one direction, in 
the plane perpendicular to the vortex at its end, is, as may be deduced easily from the formula 
in section 8, exactly half of the corresponding velocity in the neighborhood of a rectilinear 
vortex filament extending to infinity in both directions. This can also be easily seen from the 
fact that two vortex filaments, each extend- 
ing to infinity in only one direction — but 
oppositely in the two cases - - form, if com- 
bined, a single filament extending to infinity 
in both directions. The total velocity 
caused at the point P by the vortex A, 

see figure 40, is / 1 > where r= y x 2 + h 2 ; its 
vertical component is 

r. x 
w K = . • • 

47ir r 


A 



Fig. 3*. 


The vertical component due to the 
vortex B is 

I\ L — x 
w B =. \ • , 

4tit r 

where r' = \(i 2 — x) 2 4- h 2 . 

Therefore the vertical component due 
to both vortices is 


r (x /, x \ 

4 7r \r J r ' 2 / 


( 31 ) 



If we assume that the lift is uniformly 
distributed over the effective span l 2J which 
again we shall take as somewhat less than 
the actual span, then, since every element of 
the wing must be turned through the angle 

(p according to the formula 4 tan <p= y, the 

direction of the air force must be turned 
also, which means a negligible change in 
the lift, but an increase in the drag of this 
wing which must be taken into account. 

It is essential then in this calculation that we pass from a condition for a monoplane to 
one in which the wing when part of a biplane has the same lift as when considered as a mono- 
plane. The angle of attack for which this condition will arise can be estimated afterwards 
from the average of the angles <p . 


Influence of one wing of a biplane upon the other; rotation of wing profile, 
alteration of its curvature. 


10. The contribution of vortex A to the increase of the drag of the upper wing in figure 40 
is evidently 

h h 



w A 
V l 


*dx = 


A, 

h 


r Cxdx _ A, r . 

4tt Vjo r 2 - \ 4tt V R> " 


r, 

r, 
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The contribution of vortex B is, by symmetry, the same. In accordance with equation (20) , 
we can put T — and thus obtain for the increase of the drag of the upper wing 


r log r 2 /r l 

13 lj 2 2irp V 2 


(32) 


By the symbol W l2 is meant that it is the drag produced by wdng 1 upon wing 2. One can 
convince himself easily that the drag W 2l which wing 2 produces upon wing 1 has the same 
magnitude. Therefore the total increase of drag due to the fact that two monoplanes which 
produce the lifts A { and A 2 are combined to form a biplane, the two lifts remaining unchanged 
(the angles of incidence of course being changed), is 


+ W lt 


i<>k 

• l 2 q & r\ 


(33) 


in which, as always q = 0 pi 72 .) 8 

Lpon the change in the magnitude of the velocity, which in accordance with the approxi- 
mation used depends only upon the disturbance velocity v in the direction of flight, only the 
I transverse vortex of the other wing has an influence. For 

any point this influence, according to our formula, is given by 


r p 

-A* "" ' 1 

V /i\ 




r f x + l ~ x \ 

1 4tt h\r r' ) 


(34) 


Fig. 40.— Velocity at a point P due to the tip 
vortices. 


in which r and r l have the same meaning as before. The 
upper wing experiences due to the lower an increase in veloc- 
ity, the lower one experiences due to the upper a decrease in 
velocity, to which correspond, respectively, an increase or a 
decrease in lift as shown by the usual formulas. If we wish to 
keep the lifts unchanged, as required in the treatment given 
above, it is necessary to change the angles of attack correspondingly. 

The effective change in the curvature 0 of the wing profile will, for simplicity’s sake, be 

discussed here only for the medial plane of the biplane, i. e., for x — It is obtained in the 

simplest manner by differentiating the angle of inclination of the air current disturbed by the 

other wing, which is, remembering that tan <p = 


Ud 

R~dy 


^tan = 


law 

Vdy 


(35) 


Outside of the vertical plane, owing to the disturbing wing, three vortices contribute to 
the magnitude of w. A side vortex contributes, at a point at the height h and the distance y 

in front of the transverse plane, a velocity v ' perpendicular to r" of the amount ^1 — 
and therefore its share of w is 

= AWi_ i_ v\ 

1 47rr 2 V r'7 - 87rr\ r" ) 


w. 


The transverse vortex contributes 


k y 

4tt r V'r' 


in which the meaning of r" and r fn may be seen from figure 41. The total w is, accordingly, 


w = 2 w x + w 3 


r V 1 _ V _ V \ 

4tt \r 2 r 2 r" r"r" f2 J 


* The mutual action of two wings placed side by side can also be calculated from the considerations stated above, and results in a decrease of 
the drag. This decrease Is of a similar kind to that which arises in the theory of a monoplane by an increase in the aspect ratio. 

9 By change in curvature of the wiug is meant that if the flow were to be kept straight and the curvature changed, the forces on the wing would 
be changed exactly as they are oa the actual wing owing to the change in the flow.— Tr. 
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The differential of this with reference to y, for the value of y = o, is, since then r'*~r and 
r'"=h 


( _ Tl t 
Wo/ - 4jrr' 


0 + l)‘ 

the curvature sought is, then, according to equation (35), 

i = r ui i\ 

II 4jt Vr\r*^h?J 


(36) 


Calculations of tlie preceding nature were made in 1912 by my assistant, A. Betz, so as 
to compare experiments with monoplanes and biplanes and to study the influence of different 
angles of attack and different degrees of stagger of the two wings of a biplane upon each other. 
The influence upon the drag was not known to us at that time, and the calculation was carried 
out so as to obtain the changes in the lift due to -w } to v and to the curvature of the streamlines. 
In this connection the change of the lift of a monoplane when flying near the earth’s surface 
was also deduced, by calculating the influence of the “mirrored wing” exactly as was that of 
the other wing of a biplane. All that was necessary was to change some algebraic signs, 
because the mirrored wing had negative lift. The theory of these 
calculations was given by Betz in the Z. F. M. f 1914, page 253. 

The results of the theory of Betz, from a more modern 
standpoint, such as adopted here, were given in the Tech- 
nische Berichte, volume 1, page 103 et seq. There one can 
find the discussion requisite for the treatment of the most 
general case of a biplane having different spans of the two 
wings and with any stagger. In the case of great stagger 
it appears, for example, that the forward wing is in an ascend- 
ing air current caused by the rear wing; the latter is in an 
intensified descending current due to the forward wing and the vortices flowing off from it. 
Corresponding to this, if the angle of attack is unchanged, the lift of the forward wing is 

increased, and that of the rear one weakened; at the same time the ratio experiences a 

decrease for the forward wing and a marked increase for the rear one. 

For a wing in the neighborhood of the ground, owing to the influence of v there is a decrease 
of lift, and conversely there is an increase of lift due to the influence of w, provided tlie angle 

of attack is kept constant, hut as the 



Kig. 41.— Curvature of w in^-profilc at 
its middle point due to velocities 
caused by transverse and tip vortices. 


result an evident decrease in the ratio 

Lift 


Owing to 

this last it is seen why in the early days of aeronautics many machines could fly only near the 
ground and could not rise far from it. Their low-powered engines were strong enough to over- 
come the diminished drag near the ground but not that in free air. 


D. THEORY OF THE MONOPLANE. 

20. If we extend the principles, which up to this point have been applied to the influence 
of one wing upon another, to the effect upon a single wing of its own vortices, it can be said in 
advance that one would expect to find in that case effects similar to those shown in the influence 
of one wing of a biplane upon the other, i. e., the existence of lift presupposes a descending 
flow in the neighborhood of the wing, owing to which the angle of attack is made greater and 
the drag is increased, both the more so the closer to the middle the vortices flowing off at 
the ends arc, i. e., the smaller the aspect ratio is. One might propose to apply the theory pre- 
viously given for biplanes by making in the formulas of this theory the gap equal to zero. Apart 
from the fact that the formulas developed do not hold for the immediate neighborhood of the 
vortex-producing wing, but must be replaced by more accurate ones, this certainly is not the 
proper path to follow, for, in the earlier treatment, we have taken the undisturbed monoplane 
as the object with which other cases are to be compared and have asked what drag, what change 
in angle of attack, etc., arc caused by adding a second wing to this monoplane. To proceed 
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according to the same method, wc must seek for the theory of monoplanes another suitable 
object of comparison. As such, the infinitely long wing will serve. Where the discussion 
previously was about change of angle of attack, increase of drag, etc., we intend now to refer 
these to the infinitely long wing as a starting point. Since in the theoretical non viscous flow 
the infinitely long wing experiences no drag, the total drag of such a wing in such a fluid must 
be due to vortices amenable to our calculations, as the following treatment will show. In a 
viscous fluid drag will arise for both wings,* infinitely long or not, which for those angles of 
attack for which the profile is said to be “good” is, according to the results of experiment, of the 
ontyr of magnitude of the frictional resistance of a plane surface. 

The carrying out of this problem is accompanied with greater difficulties than the calcula- 
tion for a biplane as given. In order to obtain the necessary assistance for the solution of the 
problem, we shall first be obliged to improve the accuracy of our picture of the vortex system. 

The density of the lift (lift per unit length) is not constant over the whole span, but in 
general falls off gradually from a maximum at the middle nearly to zero at the ends. In ac- 
cordance with what has been proved, there corresponds to this a circulation decreasing from 
within outward. Therefore, according to the theorem that by the displacement of the closed 
curve the circulation F can change only if a corresponding quantity of vortex filaments are cut, 
wc must assume that vortex filaments proceed off from the trailing edge wherever F changes. 



Kig. 42 .— Change in shape of vortex ribbons at great distances behind the wing. 


For a portion of this edge of length dx the vortex strength is therefore to be written ^ dx, and 

hence per unit length of the 

dr 

edge is * These vortex fila- 
ments flowing off, closely side 
by side, form, taken as a whole, 
a surface-like figure, which we 
shall call a “vortex ribbon.” 
For an understanding of 
this vortex ribbon we can also 
approach the subject from an 
entirelv different side. Let us consider the flow in the immediate neighborhood of the surface of 
the wing. Since the excess in pressure below the wing and the depression above it must vanish 
as one goes beyond the side edges of the wing in any manner, there must be a fall in pressure near 
these edges, which is directed outward on the lower side of the wing and inward on the upper. 
The oncoming flow, under the action of this pressure drop, while it passes along the wing, will 
receive on the lower side an additional component outward, on the upper side, one inward, 
which does not vanish later. If we assume that at the trailing edge the flow is completely 
closed again, as is the case in non viscous flow, we will therefore have a difference in direction 
between the upper and lower flow; the upper one has a relative velocity inward with reference 
to the lower one, and this is perpendicular to the mean velocity, since on account of the Ber- 
noulli equation in the absence of a pressure difference between the two layers the numerical 
values of their velocities must be the same. This relative velocity of the two flows is exactly 
the result of the surface distribution of vortices mentioned above (as the vortex theory proves, 
a surface distribution of vortices always means a discontinuity of velocity between the regions 
lying on the two sides of the surface). The relative velocity is the greater, the greater the side- 
wise pressure drop, i. e., the greater the sidewise change in lift. The picture thus obtained 
agrees in all respects with the former one. 


21. The strengths of our vortex ribbon remain unchanged during the whole flight, yet 
the separate parts of the ribbon influence each other, and there takes place, somewhat as is 
shown in figure 42, a gradual rolling up of the ribbon, as a closer examination proves. An 
exact theoretical investigation of this phenomenon is not possible at this time; it can only be 
said that the two halves of the vortex ribbon become concentrated more and more, and that 
finally at great distances from the wing there remain a pair of vortices with rather weak cores. 
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For the practical problem, which chiefly concerns us, namely, to study the reaction of the 
vortices upon the wing, it is not necessary to know these changes going on at a great distance, 
for the parts of the vortex system nearest the wing will exercise the greatest influence. We 
shall therefore not consider the gradual transformation of the vortex ribbon, and, in order to 
make the matter quite simple, we shall make the calculation as if all the vortex filaments were 
running off behind in straight lines opposite to the direction of flight. It will be seen that, 
with this assumption, the calculations may be carried out and that they furnish a theory of 
the monoplane which is very useful and capable of giving assistance in various ways. 

If we wish to establish the method referred to with greater mathematical rigor, we can 
proceed as follows: Since the complete problem is to be developed taking into account all cir- 
cumstances, we shall limit ourselves to the case of a very small lift and shall systematically 
carry through all calculations in such a manner that only the lowest power of the lift is retained, 
all higher powers being neglected. The motion of the vortex ribbon itself is proportional to 
the total circulation, therefore also proportional to the lift; it is therefore small if the lift is 
small. If the velocities caused by the vortex ribbon arc calculated, first for the ribbon in its 
actual form, then for the ribbon simplified in the manner mentioned, the difference for the two 
distributions will be small compared with the values of the velocity, therefore small of the 
second order, i. e., small as the square of the circulation. We shall therefore neglect the differ- 
ence. Considerations of this kind are capable of deciding in every case what actions should be 
taken into account and what ones may be neglected. By our simplifications we have therefore 
made the problem linear, as a mathematician says, and by this fact we have made its solution 
possible. It must be considered a specially fortunate circumstance that, even with the greatest 
values of the lift that actually occur with the usual aspect ratios, the independent motion of the 
vortex ribbon is still fairly small, so that, in the sense of this theory, all lifts which are met in 
practice may still be regarded as small. For surfaces having large chords, as, for instance, a 
square, this no longer holds. In this case there are, in addition, other reasons which prove 
that our theory 7 is no longer sufficiently accurate. This will be shown in the next paragraph. 

It has already been mentioned that the infinitely long wing will serve as an object of com- 
parison for the theory of the monoplane. We shall formulate this now more exactly by saying: 
Every separate section of the wing of length dx shall bear the same relation to the modified flow 
due to the vortex system as does a corresponding element of an infinitely long wing to the recti- 
linear flow. The additional velocities caused by 7 the vortex system vary 7 from place to place 
and also vary 7 in the direction of the chord of the wing, so that again we have to do with an 
influence of curvature. This influence is in practice not very great and will for the sake of 
simplicity 7 be neglected. This is specially 7 allowable with wings whose chords are small in 
comparison with their spans, i. e., with those of large aspect ratio. If one wishes to express 
with mathematical exactness this simplify ing assumption, it can be said that the theory 7 of an 
actual wing of finite chord is not developed, but rather that of a “lifting line.” It is clear 
that a wing of aspect ratio 1 : 6*may 7 be approximated by a lifting line, specially 7 if one considers 
that actually the lift is concentrated for the greatest part in a region nearer the leading edge. 
It is easily seen, however, that a surface in the form of a square can be approximated only 
poorly by a lifting line. 

If we assume a straight lifting line, which lies in a plane perpendicular to the direction 
of flight, the flow due to the vortices, which according to the Biot-Sarvart law, is caused by 7 its 
own elements, will not produce any velocities at the lifting line itself except the circulation 
flow around it, which would also be present for an infinitely long lifting line having the same 
circulation as at the point observed. All disturbance velocities at a point of the lifting line, 
which are to be looked upon as deviations from the infinitely long lifting line, are due therefore 
to the vortices which run off and hence can be calculated easily by 7 an integration. 

A qualitative consideration of the distribution to be expected for the disturbance velocities 
along our lifting line shows at once that — just as was the case for a biplane— the chief thing is the 
production of a descending current of air by the vortices. If we wish to retain the lift of the 
same intensity as with the infinite wing, the angle of attack must be increased, since the descend- 
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ing air stream added to the wind due to flight causes a velocity obliquely downward. In addi- 
tion. the air force, as before, must be turned through the same angle, so that a drag results. 
The rotation will be the greater, the greater the lift and the closer to the middle of the wing 

the main production of vortices is. The 



drag must therefore increase both with 
increasing lift and with decreasing span. 

A picture of what occurs with a 
wing of finite but small chord is given in 
figure 43. There the change is shown 
of the vertical velocity component along 
a straight line parallel to the direction 
of flight through the middle of the 
wing; in the upper part of the diagram, 
for the infinitely long wing, in the lower 
part, for the finite wing. We see from 
Curve I the rising flow in front of the 
wing, its transformation into a descend- 
ing one at the wing itself and the gradual 
damping of the descending component 
due to the upward pressure drop behind 
the wing. (See sec. 15.) The corre- 
sponding curve for the finite wing is 
Curve 111. It is derived from I by 


adding to the latter the descending ve- 
locity II. We recognize the rotation of 


Kig. 4T- - Win* Inn i«K finite, but small, chord. Distribution of vertical velocity 
component along a line parallel to direction of flight. 

T. Infinitely long wing. 

IT. The downward velocity prod need by vortices flowing olT. 
in. Finite wing, sum of I and If. 


the profile as well as that of the lifting 
force, which was originally perpendicular, 

through the angle <p where tan <p — y and 


w is the velocity downward at the location of the center of pressure (i. e., at the lifting line). 
If we follow the method of Lanchester, as described in section 15, the downward velocity w can 
also be looked upon as a diminution of the ascending flow at the leading edge of the wind due 
to the absence of the sidewise prolongation of the wing, i. e., to the deviation from an infinitely 
long wing which was the basis of the treatment in section 15 Dis- 
cussions very similar to this are given in Lanchester, \olume I, 
section 117. 

It may he seen from the figure that at great distances behind 
the wing the descending velocity is 2 w, which agrees with the relation 
already mentioned that the velocities due to a straight vortex fila- 
ment extending to infinity in both directions arc twice those due to 
a filament extending to infinity in one direction only, for points in the 
plane perpendicular to this vortex passing through its end point. 

22. The mathematical processes involved in carrying out the 
theory outlined above become the most simplified if one considers as 
known the law, according to which the lift is distributed over the wing. 

We shall call this the “first problem.” The calculation is made as 
follows; The distribution of lift is the circulation expressed as a function of the abscissa x. The 

strength of the vortex filament leaving an infinitely small section dx is then • dx. This 

produces at a point x', according to what has been already explained, a vertical velocity 
downward or upward of the amount 

7 __ 1 dr dx 

dW ~iir * dx * x'-x 
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In tliis x'-x takes the place of r in section 8. If the circulation falls to zero at the ends of the 
wing, as is actually the case, then all the vortices leaving the wing are of this kind. The whole 
added velocity at the position x’, assuming that the function r(z) is everywhere continuous, is 

w = } Pf - - F- (37) 

4 ir Jo ax x — x 

We must take the so-called “chief value” of the integral, which is indeterminate at the 
point x = x', i. e., the limiting value 

/ c i ’-‘ r b \ 

1”;(J, *L.) 

must be formed, as a closer examination shows. We can do this by calculating, instead of the 
value of the velocity at the lifting line, which is determined by the preponderating influence 
of the nearest elements, the value of w for a point a little above or below the lifting line. It 
is seen that this last is not indeterminate and that by passing to a zero distance from the lifting 
line it reaches the above limit. Concerning this excursus, important in itself, the preceding 
brief remarks may be sufficient. 

After the calculation of the integral of (37), the downward velocity is known as a function 
of the abscissa x' (which we later shall call x). We then also know the inclination of the 

resultant air flow, tan <p=yi the lift dA=pYV dx’ , acting on the section fix', therefore con- 
tributes to the value of the drag 

d W= tan so ■ dA=pT(x’) ■ w{x’) dx' 

since it is inclined backward by the small angle tp. The total drag is therefore 


W- 


•/>' 


) ■ w(x') ■ dx' — 


P L f P 

O Jo 


Y (x')^d£dz' 


(38) 


For a long time it was difficult to find suitable functions to express the distribution of lift, 
from which a plausible distribution of w would be obtained by equation (30* Aftei various 
attempts it was found that a distribution of lift over the span according to a half ellipse gave 
the desired solution. According to this, if the origin of coordinates is taken at the center of the 

win s> /: ‘ / * V , dr -To* 


The “ chief value” of the integral 


/ / x V i dr — r 0 j 

r = r °VwJ ’ d^b^^b j 


r 


t dt 


is equal to ir 


and therefore the integral of equation (37) is equal to an< ^ thus is independent of x 

and constant over the whole span. Hence 


F 0 

2b 


The value of F 0 is obtained from 


r*+b/» r* bn I - 

A = p V I Ydx = p V F 0 - / 1 - 

d-Uti J-WS \ 


X 

6/2 


dx=pV F, 4 b, 


giving 




4A 


77P 


Vb 


Flence 


2A 


(39) 


7Tp Vb 2 

Since w is constant there is no need of calculating the drag by an integral, for it is simply 

2 A 2 A 2 


ii) _ ® 4 _ „ _ 

"~V A ~TpV* b 2 irqb 2 


(40) 
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llio elliptical distribution, of lift, apart from its simplicity, has obtained a special meaning 
from the fact that the diag as calculated from equation (40) proved to be the smallest drag that 

is imaginable for a monoplane having given values of 
the total lift, the span and the velocity. The proof of 
this will be given later. 

It was desirable to compare this theoretical mini- 
mum dragwith the drags actually obtained. As far back 
as 1913 this was done, but, on account of the poor quality 
of the profiles then investigated, all that was done was 
to establish that the actual drag was greater than the 
theoretical. Later (1915) it was shown, upon the in- 
vestigation of good profiles, that the theoretical drag 
corresponds very closely to the relation giving the change 
of the observed drag as a function of the lift. If we 
plot in the usual manner the theoretical drag, as given 
in formula (40) as a function of the corresponding lift, 
we obtain a parabola, which runs parallel with the 
measured “polar curve’ 7 through the entire region for 
which the profile is good. (See fig. 45.) 

This process was repeated for wings of different as- 
pect ratios, and it was proved that for one and the same 
profile the difference between the measured and the theo- 
retical drags for one and the same value of the lift coeffi- 
cient had almost identically the same value in all cases. 
This part of the drag depends, however, upon the shape of the profile, and we have therefore 
called it “profile drag/’ 7 The part of the drag obtained from theory is called “edge dra«y 5 
since it depends upon the phenomena at the edges of the wings. More justifiably the expres- 
sion “induced drag 77 is used, since in fact the phenomena with the wings are to a high degree 
analogous to the induction phenomena observed with electric conductors. 

Owing to this fact that the profile drag is independent of the aspect ratio, it became possible 
from a knowledge of the actual drag for one aspect ratio to calculate it for another. To do 
this, we pass from the formula (40) for the drag to the dimensionless lift and drag coefficients, 

by letting = and we obtain then for the coefficient of the induced drag the 

relation 



Fig. 4.’>. — f’olar diagram showing ihcoretieul drag and 
obser\ cd drag. 



( 44 ) 
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llio profile drag may then he written r w „ - <\ vi . If this drag coefficient depends only 

upon the lift coefficient, then it would he evident, since it would be the difference between the 
measured and the theoretical drag coefficients, that, for the polar curves of two different wings 
having c ai - c A , c liy 

c a 2 /y (\ 2 K 

Cwt ~ rb* niy 

and therefore 


c W2 


^W1 


7 r V b, 2 



(45) 


In a similar manner a calculation for the angle of attack may be made if we presuppose an 
elliptical distribution of lift. According to our assumption there is a (“lose connection between 
the lift of the separate elements of the wing and the “effective” angle of attack, which is the 
same as the angle of attack of an infinitely long wing. This effective angle of attack, according 
to our earlier considerations, is the angle of attack of the chord with reference to the resultant 


air current. It is therefore a --=<x- <f>. If we substitute tan o y for <f>, and introduce in equa- 
tion (39) again the lift coefficient, instead of using the lift, we obtain for the comparison of two 

wings, expressing the fact that the effective angle of attack 
.s to he the same for two equal lift coefficients, the relation 



b 2 


K 

“ 3 7 r b.y 


which leads to the transformation formula 

'■) 


, C a (F; 
a '~ 7 Ah? 




(40) 


Those formulas have been found to hold for distribu- 
tions of lift which do not deviate too much from elliptical ones, although strictly speaking they 
apply only to the latter. The fact that the type of distribution does not have a marked effect 
is based upon the consideration that both in the calculation of drag and in that of the mean 
effective angles of attack we are concerned with average results. For the calculation of the dni" 
one can also introduce the thought that no quantity varies much in the neighborhood of its 
minimum. Closer investigation of the square cornered wing has shown that, if the aspect ratio 
is not too small, the lift distribution does not deviate greatly from the elliptic type, and that the 
theoretical drag for usual aspect ratios at the most is o per cent greater than for the elliptic 
distribution. As an example of these formulas we shall take 'our figures from the book published 
by the Gottingen Institute (Krgebnisse dor Aerodynamischen Yersuehsanstalt, 1 . Licferunw, 1921). 
The first and second figures show the polar curves, and the connection between lift coefficient 
and angle of incidence for seven wings of aspect ratio" 1 1 :< to 1 : 1. The last two figures give the 
results of calculating these experimental quantities from the results for the wing having the aspect 
ratio 1 : 5. It is seen that, apart from the aspect ratios 1 1 and 1:2 practically no deviations 
are present. The fact that, the square can not be correctly deduced from the aspect ratio 1 :.j 
need not excite surprise, since the theory was developed from the concept of the lifting line, and 
a square or a wing of aspect ratio 1 :2 can scarcely be properly approximated by a lifting line. 
On the other hand it is a matter of surprise that an aspect ratio of 1 3 can be sufficiently 
approximated by the imaginary construction of a lifting line. The deviations in the case of the 
square are moreover in the direction one would expect from a lift distribution expanded over 
the chord. A quantitative theory is not available in any case at the present time. 

J3. If the lift distribution is not given, but, for. example, the downward velocity, then the 
method of treatment followed hitherto may be used, by developing the downward velocity in a 
power series and determining the constants of the series given above for the lift from the constants 

‘■The American practice is to ileline aspect ratio as the ratio oi span to chord, which would involve taking the reciprocals of the ratios given 
iu the text. Tr. 
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of this power series, by the solution of linear equations. By this the lift distribution and every- 
thing else are known. 

Another method for the solution of this “second problem” will bo obtained by the fol- 
lowing consideration: The velocities at a distance behind the wing, on account of the connection 
mentioned so often between a vortex filament extending to infinity in one direction only and 
one extending to infinity in both directions, are twice as great as those in the cross section of 
the lifting line, if we do not take into account the change in shape of the vortex ribbon. We 
therefore have here, neglecting this change in shape, an illustration of a two-dimensional fluid 
flow (uniplanar flow), for which the vertical velocity components at the point where the wing 
is reached are specified. For the simple case that the vertical velocity w is constant, as was 
found to be true for the elliptical lift distribution, tixe shape of the flow that arises has been known 
for a long time. It is given in figure 49a,. It is the same as that already considered, in another 
connection, in section 15. The picture of the streamlines show clearly the velocity disconti- 
nuity between the upper and lower sides of the vortex ribbon, indicated by the nick in the stream- 
lines, and also the vortical motion around the two extreme points of the vortex ribbon, corre- 
sponding to the ends of the wing. 

Any problems of this kind can therefore be solved by means of the methods provided by 
the potential theory for the corresponding problem of two-dimensional fluid flow. We can 
not go into these matters more closely at this time; by a later opportunity some special rela- 
tions will be discussed, however. 

A “third problem'' consists in determining the lift distribution for a definite wing having a 
given shape and given angle of attack. This problem, as may he imagined, was the first we 
proposed; its solution has taken the longest, since it leads to an integral which is awkward to 
handle. Dr. Betz in 1919 succeeded after very great efforts in solving it for the case of a square- 
cornered wing having everywhere the same profile and the same angle of attack. The wav the 
solution was obtained may be indicated briefly here. We start, as before, from the relation 


a — a? -f - <f> 


V-f 


V ) 

V 


By equation (37) w is expressed in terms of the circulation. The effective angle of attack 
a f can he expressed in terms of I\ since, according to the assumptions made before the lift, 
distribution, which is proportional to I\ depends directly upon a . The relation between a ' 
and T can he assumed to be given sufficiently exactly for our purposes by a linear expression 

r = 17 (cj ex +c 2 ) (47) 


in which is the length of the chord (measured in the' direction of flight). Bv the introduction 
of the factor 17, c x and c 2 are made pure numbers. The numerical value of c 1} which is the more 
important, can he expressed, if c ax is the lift coefficient for the infinitely long wing at the angle 
of attack a' y by the relation 

» _ 1 < 

,- 2 da’ 


In fact 


Ca , rr 
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Iuf 


It. 


P V VI 

1 T/i 

o p 1 



(/*!«' -4- <■->) 


For a flat-plate theory proves that (\ = 7r, for curved wings it has a slightly greater value. 

If, according to what has gone before, we express a bv T and ir by and write 

" ax 

^=/(x) and therefore V — £ f{x)<1x 

we obtain after a simple calculation the integral equation 

I f(x)dx + < . lt [ = 17 (Cja + Co) = const. 

Jo 4 ttJ 0 x — x 


( 48 ) 
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A solution of this equation can be obtained by expanding V as in equation (41) and develop- 
ing then all the expressions in power series of £ = ^ 9 * For every power of £ there is then 

a linear equation between the quantities r 0 , r 2 , etc. There is a system of equations, then, 
with an indefinite number of equations for an infinite number of unknowns, the solution of 
which in this form is not yet possible. The aspect ratio of the wing appears in these equations 

as a parameter, and it is clear that the solution for a small aspect ratio, i. e., ^ is easier than for 


a large one. Dr. Betz proved that a development in powers can be"made for the unknowns in 
terms of a parameter containing the aspect ratio. The calculations which are contained in the 

dissertation 11 of Dr. Betz (1919) are 



very complicated and can not be re- 
produced here ; but certain results will 
be mentioned. The Betz parameter 
L has the meaning 

2b da r 
~~c x t~ t <Zc a * 

In the application to surfaces which 
are investigated in wind tunnels the 

value ( 7 Cr is known, not ^f a r* For 
da ua 

this case theory gives a relation which 

can be expressed approximately 

3.85 b . f® — 1.3. 
t dc & 

We can thus obtain the value of -r— - 

from the connection mentioned. 

The distribution of lift density 
ever the span is elliptical for very small 


Fig. 50.— Cluing' in dislrit.ution of lid, as a fund ion of L, the parameter of Betz. 


aspect ratios and for greater ratios 
becomes more and more uniform; for 


very elongated wings it approaches gradually a rectangular distribution. Figure 50 shows 
this change in the distribution depending upon L. 

The drag of the wing with the rectangular distribution is greater naturally than with the 
elliptic distribution, since this gives the minimum of drag, yet the differences are not very great; 

for instance, for L = 4, ^ ^0^ it is about 5 per cent greater than that of the elliptical distri- 
bution. An approximation formula, according to the values obtained by Betz, is 



(0.99 + 0.015 L). 


This is applicable for values of L between 1 and 10. 

The distribution of lift, downward velocity, and drag upon a very elongated wing is shown 
qualitatively in figure 51. It is seen that the downward velocity and the drag gradually accumu- 
late at the ends of the wing. This gives also the correct transition to an infinitely long wing, 
with which for interior positions the lift is constant, and the downward velocity and the drag 
are equal to zero, while, as we know, near the ends these last quantities always assume finite 
values. 

u Printed in extracts in Heft 2 of the “ Beriehte u. Abhl. der Wiss. Ces. f. Luftf.” Munich, 1920 ( U. Oldcnbcurg). 
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E. IMPROVED THEORY OF AIRPLANES HAVING MORE THAN ONE WING. 

24. The knowledge obtained in the theory of a monoplane can be applied also to multi- 
planes and furnishes here a series of remarkable theorems. We shall limit ourselves to the 
theory of the first order, as designated in the theory of monoplanes, therefore we shall neglect 
the influences of v. Further, we shall not take into account the effect of curvature i. e* we 
shall consider the separate wings replaced by “lifting lines/ 7 For the sake of simplicity we 
shall limit ourselves to multiplanes with wings which are straight and parallel to each other. 
The generalization of the theorems for nonparallel wings, corresponding to the deduction given in 
\\ ing theory II,” will then be stated without proof. 

Let us first solve the introductory problem of 
calculating the vertical velocity w produced by a 
lifting line at a point A which lies off the lifting 
line. At the beginning let us assume that this point 
lies in the same “transverse plane 77 (plane perpen- 
dicular to the direction of flight). According to 
our assumption as to the location of A, the action 
of the transverse vortex is zero: With reference to 

the longitudinal vortices it is to be remembered that 

the velocity ^ ^ ^ produced by a longitudinal 

dl 




vortex 


of strength 


dx 


dx is perpendicular to the Distribution of lift, down-wash sincHrap for a long 


line a (see fig. o2), and therefore must be multiplied by sin /3 to obtain the vertical component. 
We arrive at the downward velocity, therefore, by integrating over the lifting line, 


viz : 


w — — 


1 I'ir ^ 

*7r Jo ax \ a / 


(49) 


This relation can be brought into another form by a partial integration. Since at both wing 
ends T = o, we have 


But 


so that we have 


d /sin 0\ d / \ a- - 2x- = 1 - 2 si 
lx \ a ) dx. \<r / a 4 ~ a 2 

| r < t>s 2 2 1 
Jo « 2 


sin 2 d 


cos 2 (3 
a 2 


4tt J 0 a a 2 u * (50) 

With the aid of this relation we can write down immediately the value of the drag which arises 
owing to a second wing being under the influence of the disturbance caused by the first wing 
tying ln the same transverse plane. Let us call vj ]2 the disturbance velocity at 
a point A on the second wing. According to the results of section 22, the 
drag then is 

fb. 

Hp = p Yji'ydx 

Jo 

as given by equation (50) is substituted, 

cos 2/3 
a 2 



or, if the value of 


Fir;, Velocity at 
a point. A off the 
lifting line, but in 
the transverse 




(51) 


The double integral, as one sees, is perfectly symmetrical in the quantities asso- 
pian" r due to the ciatcd with both wings 1 and 2. We conclude from this that the drag which 
wing 1 experiences owing to the presence of wing 2 is of the same amount as 
the drag calculated here, that, therefore, 


f I2 = w n . 
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In the more general case of two curved lifting lines lying in a transverse plane, a formula is 
obtained which differs from equation (51) only in having cos (0, + 0 2 ) in place of cos 2 0 in which 
0, and 0 2 are the angles which the line a makes with the normals on the two lifting elements con- 
nected by a, and in having ds^ls., in place of dx^dx^. The relation 11,2= ^ -m therefore holds in 
this case also. This mutual relation, which was discovered in a different manner by my assist- 
ant, Dr. Munk, is of importance in various applications. Since it plainly is not necessary for 
the lifting elements taken as a whole to belong to a single surface, the theorem may be stated : 

If, out of a lifting system all of whose elements lie in a transverse plane any two groups 
are selected, the portion of the drag experienced by group 1 due to .the velocity field of group 2 
is exactly of the same amount as that experienced by group 2 due to the velocity field of group 1 . 

We can interpret the partial integration performed above by saying that the velocity w 
appears by it as built up out of the contributions by merely infinitesimal wings having the 
length dx and the circulation T, while previously we have always built it up out of the actions 

of the separate vortices The integrand of equation (50) in fact agrees with the velocity 


which is caused by two vortex lines of equal but opposite strengths F lying at a distance dx apart. 
The double integral in equation (51) can, from this point of view, be looked upon as the sum 
of the actions of the vortex strips of all the elements dr, on all the 
lifting elements dx 2 . 

The objection might be raised that equations (50) and (51) cease 
to be applicable if the value a^o appears, since this gives an expres- 
sion of the form cc — go . They are not, therefore, suited for the 
calculation of the velocity w at the wing itself. In this ease we must 
return to equation (49), and take the “chief value’' of the integral; 
or, the value of u\ 2 and of W l2 for a lifting line that lies very close 
must be calculated, and then we can obtain our final result by passing 
to the limit for coinciding lifting lines. As is seen from this, the 
relations W l2 = lb 21 hold also for lifting linos coinciding in space, 
which, besides, may have any arbitrary lift distribution. 

The mutual drag need not, as has already been mentioned, always be positive. For in- 
stance, it is negative for two wings placed side by side, since then each wing is in an ascending 
current caused by the other, and the total drag is therefore less than the sum of the mutual 
drags which each of the wings would have at a greater distance apart. The behavior of certain 
birds which in a common flight space themselves in a regular phalanx can be explained by 
reference to this. 



25. In order to be able to treat the case of staggered wing systems, the next problem is 
to calculate the velocity field due to a lifting element of the length dx together with its pair of 
vortices at a point A which may now lie of! the transverse plane, and at a distance y from it. 
(See fig. 53.) The origin of coordinates will be taken at the projection of the point A upon 
the transverse plane, and the X axis parallel to the direction of the element. Using the abbre- 
viations 

a 2 — x 2 -f 2 l \ r 2 = a 2 + y 2 j 


the velocity produced at the point A by one of the two vortices, by formula (6b), is given by 

r 

4tt a' 

selves, is, then, putting sin 0=-* 


the component in the direction of the Z axis, to which wc here again limit our- 

U ' = w(' + r) 

The pair of vortices produces then a velocity which may be written as the difference of the 
effects of two vortices which are close together: 


dw 


dx 4tt [_( a‘ \ r / a‘r 3 JJ 
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To this must be added the contribution of the transverse vortex 

, Ydx . y 

The sum of these two velocities, if the two angles defined in figure 53 arc introduced, amounts to 

, TV/z f cos 20 , , . . sin a . cos 2 01 

dw = 4 t “T 2 <1 + sin a) i- ■ r2 j (52) 


With the help of this formula we can now calculate at once the drag experienced by a lifting 
clement situated at the point A and parallel to the former, whose length is dx, and circulation 
I' r If the first element is given the index 1, this drag is 


d 2 W.,= 


pl\\\dx l dx 2 


(1 + sin a) 


sin a cos 2 0} 


As is easily seen, the drag produced on the lifting clement 1 by the lifting element 2 is obtained 
if in place of a and 0 the values a-t-?r and 0 + ir are introducer!. J herefore it is 


„ p l\ I’, dr, dr, ftcos 2 

d W *' = 4 7r Li 


(1-sin «) - 


sin a cos* 


It is seen from this that the two parts of the drag are equal only if a = o, that is if the two ele- 
ments lie in the same transverse plane. Yet in the.general case the sum d‘ W,, + d 2 IF, is inde- 
pendent of a, therefore independent of the amount of stagger. The sum of the two mutual 
drags leads thus to the same formula as that already derived for nonstaggered wings. If we 
again pass to the general case of nonparallel lifting lines, in which again ds, and ds 2 are to be 
written in place of dx, and dx,, wc obtain as may be proved by performing the calculation, the 


relation 


wj.r 


U, ds l ds 2 cos (/Sj 


As is evident, this sum remains unchanged if the two lifting groups are displaced in the direc- 
tion of flight. Since the total drag of a lifting system is composed of such mutual drags as cal- 
culated above and of the proper drags of the separate wings, which likewise are not changed 
by a displacement of the wing in the direction of flight, the following theorem may be stated: 

The total drag of any lifting system remains unchanged if the lifting elements are displaced 
in the direction of flight without changing their lift forces. 

This “stagger theorem’’ was likewise proved by Munk. For a proper understanding of 
this theorem it must be mentioned expressly that, in the displacement of the separate lifting 
elements, their angles of attack must so be changed that the effective angles of attack and there- 
fore the lifting forces remain unaltered. . 

This theorem, which at first sight is surprising, may also be proved from considerations ot 
energy. Let us remember that, by the overcoming of the drag, work is done, and that in a non- 
viscous fluid, such as we everywhere assume, this work can not vanish. Its equivalent is, in 
fact, the kinetic energy that remains behind in the vortex motions in the rear of the lifting 
system. This energy depends only upon the character of these vortices, not upon the way in 
which they are produced. If we neglect, as we have throughout, any change in shape of the 
vortex system, then, in fact, the staggering of the separate parts of the lifting system can not 
have any influence upon the total drag. 

26. For the practical calculation of the total drag of a multiplane, we have then the follow- 
ing The total drag consists of the sum of all the separate drags and of as many mutual drags 
as” there are combinations of the wings in twos. If the nature of the lift distribution over all 
the separate wings is specified, then the proper drags are proportional to the square of the sepa- 
rate lifts; the mutual drags, to the product of the lifts of the two wings in question. If the 
coefficients of this mixed quadratic expression are all known, then one can solve without ditti- 
culty the problem: For a specified total lift, to determine the distribution of lift over the sepa- 
rate wings which will make the total drag a minimum. 
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In order to know these coefficients to a certain degree I calculated them for the case of two 
straight lifting lines whose middle points lie in the same plane of symmetry, with the assump- 
tion that the lift over each separate wing is distributed according to a half ellipse. The results 
are given in my paper “The induced drag of multiplanes ” in Volume III, part 7 of the Tech- 
nische Berichte. For this purpose the velocity w for the entire neighborhood of a wing in the 
transverse plane was first calculated by formula (49), and then the integrals for the mutual 
drags were obtained by planimetry. To show the analogy with equation (40) this may now 
be expressed by the formula 


TV, = < 


A t A 2 
7r q 6 l b 2 


(55) 


by means of which the numerical factor a can be expressed as a function of the two variables 

1 and + 2 . Calculation gave the following table: 

2 (&I+W 1 

Table I. — • Values of a 


2 ft 

&i+Aj 

0 

0. 05 

0.1 

0.15 

0.2 

0.3 

0.4 


1.0 

0.780 

0 . 655 

0.561 

0. 485 ! 

0.370 

0.290 


.8 

.690 

.600 | 

.523 

.459 

.355 

.282 

6 1 1 

.6 

.540 

.485 

.437 

.394 

.315 

.255 


The curve of the function <r is given in figure 54. For the most important case, viz, for 
two wings of equal span, I have developed an approximation formula which is 

1 — 0 . 66 ^ 

<r = \ (56) 

1.055 + 3.7 h 

It may be used from £ = 0.05 to 0.5. 

The total induced drag of a biplane is then, if b , is the greater span and if the ratio ^ 
is designated by n 

It = TF n 4-2 W n + W n = -^ 2 {A 2 4- 2o'juA 1 A 2 + v 2 A 2 2 ) (57) 

Simple calculation shows that for a given A t + A 2 this drag is a minimum for 

A 2 : A,= (m-<t) : (*-<r) (58) 

The value of the minimum is fouud to be 


tr _ (A, + A.) 2 l — a 1 

rr qb* ' 1-2?MV 


(59) 


The first factor of this formula is the drag of a monoplane having the span 6, and the lift A. + A,. 
Since <x<ju, the second factor of the formula is always less than 1, i. e., the induced drag of a 
biplane is less than that of a monoplane which in the same span carries the same load. For 
a “tandem,” i. c., an arrangement of two wings one behind the other, the stagger theorem shows 
an equivalence with two coinciding wings, i. e., a monoplane. Among the different biplanes 
having prescribed span 5, and prescribed gap h, that one is the most favorable in which the 
second wing also has the span The most favorable ratio of the two lifts is then 1 : 1 and 

the second factor of equation (59) becomes equal to * (1 +<r). 

These statements must not, howev er , be misunderstood j they refer only to the comparison 
of such wing systems as have the same value for the greatest span. Naturally, for every biplane 
a monoplane may be found with somewhat greater span than that of the biplane, which at the 
same total lift has the same induced drag as the biplane. 
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This last remark leads us to apply also to the biplane the deduction formulas obtained for 
monoplanes. All that is necessary is to replace the biplane of span l) x by a monoplane of a 
somewhat greater span kb lt which with reference to the drag — and, on the whole, with reference 
to the angle of attack — is equivalent to the biplane. If again we pass from the lift and drag 
to their coeflicicnts c a and c w , the formula connecting the drags of any two lifting systems 
1 and 2 is 


C Wl C^'2 — 


7T \ 2 (k 2 b 2 ) 2 ) 


in which, as is easily seen, the factor k f for a biplane having the most favorable distribution 
of lift, is the reciprocal of the square root of the second factor in formula (59). 

The tests of this formula with biplanes have shown that, when by giving a special shape to 
the wing the lift distribution was made elliptical, there was good agreement with the calcula- 
tions from monoplane experiments; with biplanes having the usual square-cornered wings, 
on the other hand, there was a discrepancy, which is to be attributed to the fact that the lift 
distribution on these biplanes deviates too far from an elliptic one. We can, however, retain 


the same transformation formula if the fac- 
tor k is determined empirically for every 
wing system; it is found to be somewhat 
smaller than according to the theory given 
above. The experiments on this point are 
not yet completed, so more accurate values 
can not as yet be given. The earlier Got- 
tingen experiments w T ere worked up by Dr. 
Munk, to whom this last idea is due, in the 
paper “Contribution to the aerodynamics 
of the lifting parts of airplanes” in the 
Technische Berichte, Volume II, page 187. 

27. In the previous section I have treated 
the problem of finding the minimum of the 
induced drag of a multiplane, under very 
definite assumptions concerning the distri- 
bution of lift over each separate wing. The 
strict minimum problem is however differ- 
ent, viz: 

To determine for a given front view of 



a lifting system that distribution of lift over 

all the lifting elements which will make the induced drag a minimum for a specified total lift. 

In this statement of the problem the expression “a given front view”— i. e., more exactly 
stated, a given projection of the lifting system upon a plane perpendicular to the direction of 
flight — is used to mean that the wing chord is of secondary consideration, and does not need 
to be determined until later when the selection of suitable angles of attack is made. 

The general solution of this problem was also given by L)r. Munk. It will be deduced here 
in a simpler manner than that given in Munk’s dissertation, where the solution was obtained 
by the calculus of variations. By means of the stagger theorem mentioned in section 25 the 
wing system will be referred back to the corresponding nonstaggered system. For this, as we 
showed, the relation W l2 = W n holds. We shall now introduce— with the simplifying assumption 
that all the lifting elements are parallel to each other — -a variation of the lift distribution by 
adding at any one place a lift BA and at the same time taking away an equal amount at some 
other place, so that on the whole the lift, which is prescribed, remains unchanged. We must 
now consider the change in the induced drag caused by this variation. If there is superimposed 
upon the lift distribution an additional air force 5.-1 distributed over a short portion dx, there 
arises therefrom, in addition to the drag proper of the added lift — which, however, if sufficiently 
small is of the second order — a mutual drag, because on the one hand the added lift finds itself 
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in a flow having thr downward velocity w, which is due to the lifting system, and on the other 
the lifting system is in the velocity field of the added lift. The first of these two drags, as is 

easily seen, is 5/1 y ; the other party according to our theorem, lias the same value; so the total 

drag is twice this. What condition, now, must be satisfied by the sum of the mutual drags 
caused by our twofold change of the lift distribution in order to obtain the absolute minimum 
of the induced drag? I he answer is, evidently, that we will have the minimum only if by no 
change of this kind can the drag he further diminished. The sum of the induced drags, there- 
fore, can in no case he negative; also it may not be positive, because in that case by a reversal 
of the signs of the changes which we selected we could make the sum negative. Only the value 
zero is therefore allowable. Hence, if w is the vertical velocity at point 1 and w 2 that at point 
2, we have the relation 

5/1 -f S/Cj) = o 

and, therefore, since 5/1, »- — 5/l 2 , 

w x ---= u\. 

Since this holds for all the lifting elements, we have obtained the answer. The lift distribu- 
tion which in the given wing system, for a specified total lift, causes a minimum of drag is that 
which leads to the same downward velocity at all the lifting elements. With monoplanes the 
elliptical lift distribution leads to a constant downward velocity u\ W r e recognize from this 
that the elliptical distribution in fact is that distribution of lift which causes the least drag for 
a monoplane. 

The theorem can, besides, be extended easily to the case of nonparallel lifting elements 
a transverse plane. If v: u is the velocity in the transverse plane perpendicular to the 
lifting element and € is the angle between the direction of w n and that of the given total lift, 
then, as may he shown without difficulty, ■w n = u' 0 cose for all the elements. (If e — o, and 
hence cos e 1, the statement made above again appears.) 

28. A way to solve the problem of finding the lift distribution for a prescribed distribution 
of the vertical velocity has been indicated already in section 23. The velocity field left behind 
in the air by the lifting surface is, approximately, according to the remark made before, a 
uniplanar How around the vortex system produced by the lifting system in its motion, and this 
vortex system may be regarded, as a first approximation, as a solid body in the fluid. In the 
minimum case this iigure, according to the results of section 27, moves like a rigid body, not 
alone in the case of parallel lifting elements, but also in the general case, for the general mini- 
mum condition, w n = w 0 cos e, expresses directly that the normal velocity of the fluid at an 
element of the rigid figure moving in the direction of the lift coincides with the normal com- 
ponent of the velocity w 0 of the rigid figure itself. The problem is thus reduced to a perfectly 
definite one treated in the hydrodynamics of uniplanar fluid motion. 

This uniplanar flow can be brought into relation, in a specially (dear manner, with the 
pressure distribution existing on the wing system. The wing system, during its motion along 
its path, imparts to one portion of the air after the other the velocities which we have learned 
to know as the result of the vortices flowing of! from the wings. This transmission of velocity 
is the result of the spreading out of the pressure field of the wing system over the air particles 
one after another. In order to simplify the phenomenon for ourselves we can now imagine 
that these velocities are produced at the same moment by a sort of impulse phenomenon over 
the whole path of the lifting sj^stem. To produce this impulse it is necessary to have a solid 
figure of the shape of the geometrical region passed over by the lifting system (i. e., of the 
shape of the vortex surfaces which it leaves behind). If we are concerned with a system of 
least drag, this figure moves as a rigid body; otherwise it would also experience a change of 
shape due to the impulse. The final velocity w* of the figure coincides with the motion of the 
vortex surfaces at a great distance from the lifting system, and is therefore to he put equal 
to 2 w. For a monoplane having elliptical distribution our figure is therefore an infinitely 
long flat plate of the breadth b. 
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By the production of the velocity during the impulselike acceleration an increase in the 
pressure p, arises below the plate (in the case of multiplanes, under each of the plates corre- 
sponding to the separate wings) and at the same time a decrease in pressure p 2 above the plate 
(or the plates). We can now compare in a very simple manner the total action of the pressure 
differences at each point of a plate during the time of the impulse with the total action of the 
pressure differences of the wing in its forward movement at the point of the medium in question. 
If the resulting motion is the same in both cases, then the pressure differences integrated through 
the proper times must have the same values. If in the impulse phenomenon lasting a timer a 


portion of the fluid of length l is considered, and if therefore the action of the lift 


dA * 

, ax in the 
dx 


time t= y required to pass over the length l is to be compared, the following relation must 


hold for the conditions on a strip of width dx: 



lh)dt 


dA 

dx 


dx 


l 

V 


(61) 


A formula connected with our previous relations can be obtained by a tn nsformation of the 
left-hand term. According to a known extension of the Bernouilli equation for accelerated 
motion we have 


d<t> pv 2 
p dt r 2 




For our impulse phenomenon the arbitrary time function /(f) is a constant, since at the 
points of the fluid lying far away from the impinging plate the pressure does not change. If 
the impulse is sufficiently quick, then during the short time of impulse r the acceleration and 

the pressure differences will be very large, and therefore the term may be neglected in 

comparison with the other two, since it itself does not exceed moderate values. We obtain 
therefore the simplified relation 

p -b p — const. = p Q 


which, if at the beginning everything is at rest, (*!>„ = o) may be integrated to 

p<t> - j (p 0 - p)dt (62) 

we can therefore write, in equation (61), the expression p(<I> 2 — 4^) in place of The 

potential differences <t> 2 — which here appears is, according to the connection between poten- 
tial and circulation (see sec. 5), nothing but the circulation V for a closed curve which passed 
around one edge of the vortex ribbon and intersects our vortex ribbon at x, the point consid- 
ered. This circulation is again nothing but the circulation around the wing at the point x. 
If the factor Idx is omitted from both sides of equation (61), it takes, as a result of this trans- 
formation, the form 

d f x =p{%-^)V^pVY (63) 

We have thus proved in an entirely independent way, as we see, the Kutta-Joukowski 
theorem for a wing element, which previously we took over, without proof, from the infinitely 
long wing. 

The relations deduced in the previous paragraphs permit, in the case of a constant w f 
the formation of general theorems for w and W in place of (39) and (40). By integration 
of (63) the total lift is at once obtained 

A = pVzj(<t> 2 -<l’ i )dx (64) 

The values of 4> in this formula are proportional to the velocity w* of the vortex ribbon, 
that is, are dependent upon A. Quantities which are independent of A are derived if the poten- 
tials <£ are divided by w*. In this way we obtain the potentials for a velocity of the vortex 
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ribbon equal to 1. The potential, being the line integral of the velocity, has the dimension 
velocity times length; the potential <£ for «’*=1 has therefore the dimension of a length, and 
hence 

(<f> 2 — 4>\)dx 

is a surface, which will in what follows be called F' , which depends only upon the geometrical 
properties of the projection of the wing system upon a plant 1 perpendicular to the direction 
of flight, therefore upon the front view of the wing system; and which evidently for geometrically 
similar front views is proportional to the square of the span. By introducing F' into equation 
(64) we have, since = 

A= P Vw*F' (65) 


From this we may immediately deduce w* t and thereby also the downward velocity at the point 
of the -wing system 


w — 


If this value is introduced into the relation 


1 *_ A 
2 W ~2 P VF' 

T,T W , 

\\ — | rA , we 


have 


( 66 ) 




A 2 A 2 

2p V 2 F'~ 4qF' 


(67) 


The evaluation in the manner of the flow of figure 49 gives a potential <t>, if the span of 
the wing is set equal to 5, which has the value V (6/2) 2 — x 2 at the plate. The geometrical expres- 

7t6 2 

sion of this value gives a circle having the span b as diameter, therefore F f = Using this 


value formula (67) passes over in fact into formula (40). 

It may also be noted that a uniform velocity can be superimposed upon the uniplanar 
flow here discussed, whose discontinuity in potential at the rigid figure representing the vortex 
ribbon causes the surface F' , without thereby changing the relation for F', for the potential 
discontinuity between tlie lower and upper sides, with which we are here concerned, is not 
changed by the superimposed uniform motion. We may now choose the velocity of the uni- 
form motion exactly opposite and equal to the velocity w* of the rigid figure, and thereby 
secure the condition that in the new flow the rigid figure is at rest and is surrounded by a flow 
which at infinity has the velocity w. The forces whicli the rigid figure experiences by the 
production of this motion, and which are connected intimately with the so-called u apparent 
mass,” are what we have here set in parallel with the wing-lift. 

The surface F* supplies in addition a very simple mechanical connection between the 
velocity w on the one hand and the lift and drag on the other. According to equation (65) 

A = pF f Vw* 

WV= Aw = pF' r' -f 


where in the second equation use has again been made of the relation w* = 2tc. Now pF’V 
is the mass erf air flowing per second through the section F'. If in order to simplify the whole 
problem it is once assumed that all the air particles within the section F f experience the full 
deviation «?* but that all outside are entirely undeviated, then exactly the correct lift and 
the correct work due to drag arc obtained by application of the impulse theorem and the energy 
theorem. For the lift is now equal to the mass of the fluid deviated per second times the 
vertical velocity imparted to it, therefore equal to the impulse imparted to the medium. Also, 
in the same manner, the work done per second by the drag WV is the product of the mass of 
the fluid passing per second times the half square of the deflection velocity, and therefore equal 
to the kinetic energy left behind in the medium. This relation is indeed best suited to establish 
the phenomena of the theory of airplanes in a course for students who are only slightly skilled 
in mathematics. The fact that for a monoplane the circle having a diameter equal to the span 
comes out as the surface F f is one that will appear most plausible to the laity. 
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29. Of the theory pictured in the preceding section, according to which the determination 
of the induced drag in the case of the most favorable lift distribution is reduced to a problem 
of the potential theory, manifold applications have already been made. Especially, Dr. Gram- 
mel and K. Pohlhausen have treated, at my instigation, the case of the biplane made up of two 
straight monoplanes of the same span, and also, on the other hand, that of a monoplane having a 
longitudinal slot. The calculations in both cases are solved by means of elliptic integrals. I 
have given the formulas in my Wing Theory II. It may be sufficient here to state the practical 
final result, which is referred to the magnitude of the surfaces F'. These surfaces are best 
expressed for biplanes in terms of the corresponding surfaces of the monoplane having the same 

span. In fact, the ratio F f : ^ 6 2 , as is easily seen, equals the square of the factor k, introduced 

in section 26, by which the span must be increased in order to have a monoplane of the same 
induced drag. The values of Jc 2 for the biplane are obtained from the following table. The gap 
of the biplano, i. e., the distance apart of the two wings, is designated by h. 


Table 2. 

Values of fc* « F f : * ft 2 


0 . 

0.05 

0.10 

! 0.15 

0.2 

0.3 

0.4 ' 

i 0.5 

1 

1 co 1 

1.000 

1.156 

1.212 

1 . 289 

1.352 

1.461 

1.550 

1. 626 

; 2.000 


The values given in the table may be expressed by the approximation formula 

^1.027 + 3.84 h/b 
k 1 + 1.6371/6 ‘ 


( 68 ) 


In the case of the monoplane having a slot a suitable comparison wing is obtained by shov- 
ing the two halves of the monoplane together until the slot is closed. If b is the original span 
and d is the width of the slot, this monoplane has evidently the span d — d. We shall therefore 

form the ratio F' : ~ (b-d) 2 and again designate it by k 2 . Calculations gave the following 
values: 

Table 3. 


Values of k J = F': 


0.000 

0.001 

i 0.010 

i 0.0316 

0.100 

0.250 

0.500 

1.000 

1.000 

0 . 762 

! 0.676 

i 0.6200 ; 

0.568 

0.528 

0.506 , 

0 . 500 


It is seen that even very narrow slots produce an important increase in the induced drag. 
For a very wide opening 1c 2 falls to one-half, as may be deduced easily from the fact that now, 
instead of one monoplane, we have two monoplanes of half the span. The values given in the 
table may be expressed by the approximation formula 

P = 1 “ 2 Vl + 0.35 (log, 0 bid)* (69) 

Figures 55-57 show, at the left, the uniplanar w* — flow and, at the right, the surfaces F f , for a 
monoplane, a biplane, and a monoplane with a slot. 

F. AEROFOILS IN A TUBE OR IN A FREE JET. 

30. To draw conclusions from the experimental results obtained in a tube bounded by 
solid walls or in a free jet from a nozzle, it is very useful to know the influence of the neighbor- 
ing walls and of the boundaries of the jet upon the phenomena at the aerofoil. We wish indeed 
to know the behavior of the aerofoil in an air space infinitely extended in all directions; and the 
problem therefore arises to introduce a method for passing by calculation from the case which 
prevails in the experiments to that of the unlimited air space. For this purpose we shall next 
state clearly the boundary conditions which exist at solid walls parallel to the direction of the 
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wind and at the free boundary of a jet. At solid walls, the velocity components normal to the 
wall,u* n , must equal zero; on a free jet boundary, on the other hand, the pressure is to be put 
equal to that of the surrounding quiescent air layer, and therefore is constant. We can trans- 
form these last relations as follows, keeping within our theory of the first order. According to 
the Bernouilli theorem, if l 7 is the undisturbed wind velocity, and u, v, and w are the additional 
velocities 


or, since p — p li: 


p + 2 [A + ( V+ v) 3 + w 2 J = p 0 t- g V 3 
u 2 + v 2 4- w 2 -1- 2 Vv — o. 


If we neglect the squares of the disturbance velocities as being small of the second order, we 
have as the approximate initial condition for the free jet v — o. We proceed a step farther upon 
the path indicated to us by the approximation theory of the first order if we prescribe the value 
v = o, not for the actual jet boundary, but for that cylinder which is given by the' surface of the 
undeviated jet. By doing this the boundary condition for the free jet becomes very similar, 
in a formal way, to that for the solid walls. 

The two problems can now be solved in the following manner: We consider first the velocity 
field for the unlimited air space, according to the exposition previously given. This field offers, 
both in the case of the tube and in that of the free jet, contradictions with our boundary condi- 
tions at the walls or the jet boundary. We must superimpose a velocity field which in the 
interior of the region considered is free of singularities and which on the boundaries has veloci- 
ties opposite to those velocity components, the vanishing of which is prescribed by the boundary 
condition. It is easily seen that by the superposition of this second velocity field on the original 
one the boundary conditions are satisfied exactly. The influence of this second field upon 
the aerofoil is now exactly that influence which we are seeking, and which we can calculate 
from the results of the theory of aerofoils as soon as this second field is known. 

The additional velocity field corresponds to a pure potential motion; we have, therefore, the 
problem of determining its potential In the case of solid walls we are thus led to the problem 
of finding the potential for a given region (the interior of the tube) when the normal component 
v) n of the flow is given at the boundary of the region. This is the so-called “second boundary 
value problem” of the potential theory. The corresponding problem for a jet, as we shall see 
at once, leads to the “first boundary value problem,” in which at the boundary the values of 
the potential itself are prescribed. According to what has been said above our region is a cylinder 
whose generating lines are parallel to the velocity V, hence, parallel to the Y axis, and for each 


point on the boundary the relation ^ = — v (in which the dashes indicate boundary values) is 
prescribed. Integrating this relation for each generating line gives 


(?/) = ~ P vdy 

J — no 


If we go sufficiently far upstream every influence of the aerofoil vanishes; therefore for y -= — oc, 
$ = o; and hence y= - go is taken as the lower limit of the integral. By this, then, we obtain 
the boundary values of the potential $ (?/). 

The complete calculation of the added potential for the entire interior of the tube or jet is 
fairly difficult. If we concern ourselves, however, only with our main problem, to determine 
the corrections which must be applied to our experimental results, then we can again assume 
that the velocity components perpendicular to the axis of the tube in the plane of our aerofoil 
are half as large as at a great distance behind it. This consideration, which proceeded from the 
comparison of a vortex filament proceeding to infinity in one direction only with one proceeding 
to infinity in both directions, holds here exactly as well as in the cases discussed previously. 
We can therefore pass here as before from the space problem to a uniplanar one if we calculate 
the phenomena far behind the aerofoil. Our boundary conditions for the uniplanar problem are, 
for the tube, w n - e, for the free jet, 3> = const. The last condition may be interpreted specially 
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conveniently if use is made of the method of treatment of section 28. Since in this, for the end 
of the impulse phenomenon, $= I {p 0 — p) dt, <t> — const, means simply that p = const., which? 

1 Jo 

indeed, was the original boundary condition for the free jet. 

31. The conditions stated in the preceding section can be secured most easily for a jet, or 
tube, of a circular cross section. In this case the added motion is obtained very simph r by 
assuming for every vortex flowing off an equally strong one outside the circle, at the point out- 
side corresponding to the one inside according to the reciprocal radii. If the direction of rota- 
tion of the external vortex is taken the same as that of the interior one, then at the points of 
the circle the boundary condition for a free jet is obtained; and, if opposite directions of rotations 
are taken, then the boundary condition for a tube is satisfied. This may be expressed by saying 
that there is combined with the aerofoil another obtained by reflexion according to reciprocal 
radii, whose circulation at corresponding points is the same in absolute value as that of the 
actual aerofoil, anti for the jet it has the same sign, but for the tube the opposite sign. 

The exact calculation has been made for a straight monoplane in the middle of the jet, 
assuming the lift to be distributed according to a half ellipse. If b is the span of the monoplane 
and I) the diameter of the jet, then the disturbance velocity w r caused by the jet boundary at 
the distance x from the middle of the jet is 


w — 


>+;*■ 


is'+^r+ote 


w= 


A ( 

7 rl) 2 p \ \ 

>m this 

42 l“i + V*Y+ 

3 p PL 16 \d) h 64 \DJ 


•) 


(70) 


in which £ = 2 xb/l) 2 . 

The added drag calculated from this velocity according to equation (38) is found to be 

A : 


Ti ) + ■ 


] 


(71) 


A similar calculation for a uniform lift distribution gave for the first term in the formula 
for the drag the same value as in equation (71). It appears that the other terms of the series 
have but little importance with the usual ratios, so that we can limit ourselves to the first term. 
An approximation treatment shows, further, that any small wing system, in the middle of the 
circular jet gives rise to the same expression. We can therefore write for the total induced 
drag of the wing system in a jet of cross section F 0 , if the surface F T is again introduced from 
section 28, 

A 2 


W= 


4 q 


G" + 2 V 0 ) 


(72) 


For a tube of circular cross section the same disturbance effect is found, but with the oppo- 
site sign; and therefore we have the approximation formula for the drag 

A 2 


W=-( 1 - 1 ^ 
4q\F f 2 FJ 


(72a) 


The correction, owing to the consideration of the finite cross section of the jet, is for the 

ratios ordinarily used not small. For ^ = \ > it is already one-eighth of the induced drag. 

Formula (71) gives 0.1262 instead of 0.125; the corresponding formula for uniform distribution 
gives 0.127. It is seen, therefore, that the differences are not great, and that the approxima- 
tion formula (72) is satisfactory for most cases. 

For a tube of rectangular cross section the calculations would have to he made in such a 
manner that the aerofoil was mirrored at all the walls an infinite number of times, like a check- 
erboard. Further development of the calculation leads to elliptic functions. It has not yet 
been carried through. One can assume, however, that for a tube having a square cross section 
the influence of the walls will be of about the same magnitude as for the circle having an equal 
area. 
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G. APPLICATION OF THE THEORY OF AEROFOILS TO THE SCREW PROPELLER. 

32. The fundamental ideas of the aerofoil theory can be applied step by step to the screw 
propeller. For the elements of the blades the Kutta-Joukowski formula holds, viz, that the air 
force is perpendicular to the velocity c of the element with reference to the air and 'that, per 
unit length of the blade, it has the value pTc. Corresponding to what has gone before, vorti- 
ces will arise at the blade, having a vortex strength per unit length equal to If we wish 

again to construct a theory of the first order, that is, if we agree to consider as small the air 
forces and the velocities produced by them, then again the proper motion of the vortices will be 
small and therefore in a first approximation may again be neglected. The vortices then have 
the shape of screw lines and form vortex ribbons which — if for the sake of simplicity we assume 
straight radial blades — have the shape of ordinary screw surfaces. 

The calculation of the velocity field of a screw vortex is markedly more complicated than 
that of a rectilinear vortex and leads to functions which thus far have not been studied in detail. 
In spite of this it is possible, as Dr. Betz has shown, to prove a series of general theorems very 
similar to those of Munk for multiplanes. Since the velocity c is not the same at the separate 
blade elements, we must speak of the u work lost ” where Munk speaks of drag. The work applied 
for the motion of the propeller is composed of two parts — useful work -f work lost. The latter 
in our ideal case, where friction is excluded, is transformed completely into kinetic energy of the 
air. The kinetic energy stands again in close connection with the vortex system produced by 
the propeller. Betz proved, among others, the following theorems: 

(1) If two elements of a propeller blade lie upon the same radius at distances x and £ from 
the axis, then the work lost at the point £ due to the disturbance velocity caused by the air 
force at the point x is equal to the work lost at the point x owing to the disturbance velocity 
caused by the air force at the point £. 

(2) This theorem must be somewhat modified for two elements which do not lie on the same 
radius. It reads: The work lost at the point £ due to the disturbance velocity caused by the 
air force at the point x is of the same amount as the work which would be lost at the point x if 
the screw vortex proceeding from the element at £ were to pass out forward in the prolongation 
of the actual vortex instead of going backward. 12 

(3) This last theorem leads at once to the following relation for the sum of the two amounts 
of work lost: The total work lost due to the mutual action of the airforces by two blade elements 
at points x and £ is the same as the work which would be lost at one point alone if the screw vor- 
tex proceeding from the other point were to extend to infinity both forward and backward. 

It is easily seen that this theorem is perfect^ analogous to the stagger theorem of section 
25, for if the vortex of the inducing element extends in both directions, then the position of the 
element itself on its own vortex strip is immaterial as far as the velocity field produced is con- 
cerned. 13 It is therefore true of screws that nothing is changed in the total energy-loss if blade 
elements are displaced in any way, without change of their air forces, along the relative stream- 
lines passing through them (i. e., in this case, screw lines). This naturally is connected again 
with the fact that the total amount of the energy loss depends only upon the final distribution 
of the vortex systems, not upon the relative position of the places where the separate vortices 
arise. 

Theorem No. 3 will be of use to us also in what follows. It can be made clearer by the fol- 
lowing consideration. The field of the vortex ribbon of a lifting element dies away very quickly 
forward of the element, but in the rear it extends over the entire length of the path traversed. 

If the sum is formed of the two mutual losses in work of two elements at the points x and£,we 
can proceed, owing to the stagger theorem, to displace one of the two elements along its screw 
line so far backward that its velocity field is no longer appreciable at the position of the other 

« fn this the sense of rotation of the transverse vortex is to be reversed. 

* a The transverse vortex in this case cancels out completely in the determination of the velocity Celd, since it appears twice with opposite 
senses of circulation. 

20167—23 14 
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undisplaced clement. At the same time, however, the influence of the hitter element upon the 
iiv<t is increased since its vortex ribbon, viewed from the new position of the first element, extends 
as far forward as backward. The sum of the two mutual losses in work is reduced in this manner 
to the loss which the velocity field of the front element produces upon the displaced one. 14 

(4) The most important of Betz’s theorems, from a practical standpoint, furnishes the com- 
plete analogy to Munk’s theorem concerning the wing system having the least drag, and, corre- 
sponding perfectly to the statements in sections 27 and 28, maybe expressed thus: The flow 
behind a propeller having the least loss in energy is as if the screw surfaces passed over by the 
propeller blades were solidified into a solid figure and this were displaced backward in the 
1 nonviscous fluid with a given small velocity. The potential difler- 

encc between the front and rear sides of a screw surface at one and 
\ the same point furnishes, then, again the circulation V of the corro- 

\ A sponding point of the propeller blade. 

\ r \ A short proof of theorem 4 will be given. For this purpose the 

principal equations for the action of a screw must first be deduced. 
Tl u » -screw is imagined to be displaced with the velocity v relative 
w to the air, and to rotate at the same time with the angular velocity a?. 

[ jg. rs. — velocity field near :i bi-.vic \ blade element at the distance x from the axis has then, with ref- 

ofascrevv-proreiier. erencc to the air which in the theory of the first order may be 

assumed to be at rest, the velocity c 7 with the components v and xx. (See fig. 58.) 

If no vortex were produced, then, with the assumption of a nonviscous fluid, an air force 
dP would arise, which, according to the Kutta-Joukowski theorems, would be perpendicular 
to the velocity c and would have the value, for a blade element of length dx f 

dP = pVcdx (75) 

The force dP is decomposed into two components, of which the one in the direction of v 
interests us specially, since it is applied to the screw. This component is 

dS = dP . cos e pTxxdx (74) 

The total thrust, if there are n blades, is then 


The other component 


dT = dP sin € - pYv dx 


furnishes a contribution as a torque to the rotation moment. It is seen at once that dS . v = 
dJ x «, i. e., the useful thrustwork is equal to the work done by the torque hitherto used 
in our calculations. This depends immediately upon our assumption that the force dP is 
perpendicular to the velocity c. But the screw blades actually produce a vortex system and 
we must ask as to the reaction of the vortex system upon the phenomena of a screw. We 
shall assume, exactly as in the aerofoil theory, that we turn the blade profile in such a manner 
that the lifting forces desired by us actually come into play. Since we are interested here 
merely in the loss in work caused by the vortex system, we have to do only with the drag 
components caused by the vortex system. This depends, exactly as before, upon the velocity 
component perpendicular to the velocity of motion of the element, wliieh in tills case equals c. 
We shall again designate it by w. The added velocity component w furnishes a drag in the 
direction of motion equal to 

(IQ — dP ™ =pVw tlx (77) 

The loss of work per second is therefore 

dQ . c = pYdx . w . c ( = d-P . w) 


u Tllis process of thought can be applied, naturally, in the same way to aerofoils ami furnishes a convenient deduction for the sum of the 
drags Wi«+ W n . 
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Tf now, according to figure 58, we put c v/s\n E, our problem is to make a minimum the total 
loss of work 


n f I ' w< 
sin 


The variation of this quantity must therefore be put equal to zero. We shall proceed, for 
this purpose, similarly to the way Munk’s theorem was deduced in section 27. We shall 
change by small amounts the circulation at two places, which may reach from to x 
and fro mV to in such a manner that the total thrust remains unchanged. According 

to equation (75) we must make 

5\\ . jy/jvl- 51', . xJx. r =o 

Kxactly as before the condition for the minimum is obtained if the loss of energy due to our 
added circulation remains unchanged. In order to calculate the loss, let us make use of 
theorem No. 3 and assume that the added wing forces are brought into action far behind the 
propeller so that the loss is merely the product of the added air force by the velocity w* 
which arises from the vortices of the propeller, and therefore for the first element is equal to 


pvdV. - 1 dx r 
y 1 sin e 1 


Omit tin" the constant factor, we obtain as the minimum condition 


in * * 

5l\ 1 dx^dV, • 2 

1 sin e, 1 1 sin 


from which is derived, making use of equation (79) 


= const. 


aq sin x 2 sin e 2 

We must compare this condition with that obtained for the velocity components normal 
to a rigid screw surface, when this surface is moved backward with the ^ r ^ 

velocity w\ We then have (sec iig. 59): , 


But on a screw surface the pitch h is connected with the angle of pitch 
and the radius x by the relation h = 2ttx tan e. 

Multiplying this last equation with that for ?r n , and solving, we get 

2t tw'x sin e 
Wn= h 


and therefore 


= const. 


On comparing (81) with (80) it is seen that by a suitable choice of w' the 
value of w a can always he made to agree with that oi w*, which proves 
Bet// theorem. 


Fio. 59. — Ide-ili/ecl vortex sys- 
tem of a screw-prorellpr. 


33. In order to learn more accurately the nature of the distribution of circulation which 
we are seeking, we shall proceed as if the velocity field at great distances from the screw is pro- 
duced by having the velocity w f in the direction of the axis imparted impulsively to the rigid 
liguro composed of the screw surfaces. In a purely qualitative way out* can sec that with any 
system of screw surfaces having a small pilch the air in the interior of the system is actually 
accelerated backward, with, of course, the appearance of tangential velocity components whose 
intensity is a function of the angle e and is greatest for t-=45°. At the axis itself there is neither 
an axial nor a tangential acceleration. Less simple are the conditions near the outer boundary 
of the screw surface where a How around the edges of the surfaces occurs. 

In order to obtain a quantitative statement, we shall for simplicity’s sake next think of a 
screw having a large number of blades. Our rigid figure consists, then, of a very large number 
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of screw surfaces, lying dose together, and therefore the air is led with difficulty into the interior. 
When the impulse occurs it escapes in the direction of the normals to the screw surface. The 
radial velocity components iv r will be appreciable only in the neighborhood of the outer boun- 
dary of the screw surfaces; further in, we may put it equal to zero approximately. For the 
tangential components w t and the axial components w- a , the relations hold, as is easily seen, 

w t ~w* sin t — w* cos e sin e 
iu a = W* COS € = w' COS 2 e. 

The angle may be expressed by writing 


f h 

tan e = -=■ — 
2irx 


V 

Xu 


In this, for brevity’s 
screw, tan «, =1). Then 

and hence 


(82) 

Ji V 

sake, 0 = is put equal to r' ( r ' is that radius for which the pitch of the 


sin e = 


T X 

, , and cos e = 

(83) 

\r 2 -j-x* 2 + T" 

>' • , 2 ? and w H = w' , 7 2 

r' 2 +x 2 H r' 2 +x 2 

(84) 


, (='?) 


(85) 


We must now determine the circulation around the separate blades as a function of the 
radius x. For a screw with n blades the total circulation of the vortices inside the circle of 
radius x coincides with the line integral for the closed circle of radius x ; this circulation must 
evidently equal nY, where r is the circulation of one of the screw blades at 
the point x. From this we have 

2t tx • w t = 2tt r'w' x 2 

n n r' 2 -fr 

The curve for F according to equation (So) is shown in curve l of figure 62. 

At the ends of the blades wc would expect to have a decrease of the cir- 
culation of the same character as found for aerofoils. An approximate treat- 
ment can be devised in the following way: We imagine an infinite series of 
aerofoils which have a distance apart a and are not staggered and which 
extend infinitely far toward the left. We inquire what is the most favor- 
able distribution of lift near the ends of these aerofoils. The distance a is 
then to be made equal to the perpendicular distance apart of the edges of 
two consecutive blades of the screw, i. e., according to Figure 59, 



Fki. t»0. — Conformal trans- 
form at ion 


•SO*})- 


h 

a = cos e, 
n 1 


27T7’ / 

n 


r - 

r ' 2 r 2 


( 86 ) 


The problem now, according to the procedure of section 28, may be solved by seeking the po- 
tential flow around the edges of the corresponding family of planes and by determining the dis- 
continuity of potential at the planes. This problem may be solved without difficulty by means 
of conformal representations. (See see. 10.) It can be shown that the plane with the straight 
cuts as shown in figure 00, which we shall call the z plane, may be transformed into the unit 
circle (t plane) by the formula 

a \ / i \ 

(87) 


a 1 

2 = lo< 


2 ( < + 0 


The flow of figure 58 is transformed thereby into circulation flow around the unit circle; in fact 

<i> +i ^ = iC log t 

After a short calculation, by elimination of t, we have 


a . + / T 

z=- log cos ^ 

7 r B C 


(88) 
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For the surface around which the flow takes place, and which is given by the streamline \ p = o, 
we have, therefore, 

a , <I> 

X= log COS y, 

7T O 

or 

4> = ± C cos -1 e a (89) 


which gives the real values for negative values of x. The potentials thus obtained or the 
velocity equal to 1 of the free flow ^to obtain which C must be put equal to which, accord- 
ing to what has gone before, give us the surface F' } form a picture such as is shown in figure 
61. Bv means of this one can form a definite judgment as to how the circulation, and with it 
the thrust also, decreases at the blade tips. We can replace the shaded portions of Figure 61 
by a straight line, having an equal area below it which, in accordance with 
the integration performed, must lie behind the blade tips at the distance 

a' = a log 2 = 0.2207 a (90) 

We conclude from this that, with screws also, the decrease of circulation 
at the blade tips has about the same effect as if the screw had a radius dimin- 
ished by 0.2207a, and then the air would be considered uniform in every circle of 
radius x (as would be the case for a screw having an infinite number of blades). 

The properties found for the inner portion of the screw and for its edge 
may be combined into a single formula, which can be applied as an approxima- 
tion formula also for screws having a small number of blades. This formula is obtained by 

2 -*■ O’— x) 

multiplying the value in equation (85) by the expression “ cos -1 e a } which for large values 



Fig. 61. — Potential ol> 
tained by the flow of 
figure 60. 


T — X 

of takes the value 1. 
a 

Thus we obtain the formula 


(For — x of formula (89) r — x is here substituted, as is obvious.) 


_ 4 w/r' 
n 


x 2 

r' 2 j rX 2 


-7T (r-x) 

cos -1 e a 


(91) 


The curve of F according to equation (91) is given in figure 62, for a 4-blade screw and for 
r': r — 1 : 5, which correspond to average conditions in practice. 

The whole deduction holds, as has already been remarked, for screws which are not heavily 
loaded. For screws with heavy loads an improvement can be intro- 
duced by calculating the pitch of the screw surfaces formed by the 
vortices, corresponding to the state of flow prevailing in the circular 

v 

plane of the screw. Instead of writing tan £ = , we must write, more 


2 

exactly, tan £ — , in which V is the velocity of flight, since in 

Wt 
2 

(]io screw disk piano half of the final disturbance velocities is already 

iv f 

present. A useful approximation is obtained if, retaining our formulae, v is put equal to 1 f 2 > 

w' 

and therefore r' is put equal to F+ 2 . 

ui 

After the circulation is known, the distribution of thrust and torque may be calculated 
easily by means of equations (74) and (76), and thus, following the method used in the aerofoil 
theory, the requisite widths of the blades and angles of attack may be determined in order that 
for a given working condition (i. e., r' and w' given), in which the screw is to have the most 
favorable performance, all the information may be deduced from the theory. By taking into 
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Fh;. 62.— Distribution of circu- 
lation. I. Tnfinilo number of 
blades. II. Four blades a is 
the distance of two vortices. 
(See Fig. W.) 
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account the more exact velocity relations in the propeller-disk plane this information may 
be improved. 


The aerofoil theory has numerous further applications. An investigation of curved flights 
specially of the moment— important in discussions of stability— around -the longitudinal axis in the 
case of a wing moved in a circle, is at present being made, also the calculation of the moment 
of a warped wing. A series of not unimportant single questions must wait for a further im- 
provement of the theory, c. g., various conclusions specially concerning properties of profiles, 
influence of curvature, etc., can be reached, if we pass from the lifting line to the case of a load 
distributed also along the chord; for the treatment of a wing set oblique to the direction of 
flight the assumption of a load distributed along the chord is necessary since in this case the 
conditions contradict the “lifting line.” Investigations of this kind, which can be accom- 
plished only by very comprehensive numerical calculations, were begun during the war but 
since then, owing to a lack of fellow workers, have had to remain unfinished. A similar state- 
ment also applies to 'the calculations of a flapping wing already begun, in which one is likewise 
forced to assume the lift distributed along the chord, since otherwise the result is indefinite. 
Therefore much remains to be done. 


MOST IMPORTANT SYMBOLS. 

p =density. 

V = velocity of ihe airplane. 

u, v, a; = velocity components in the X, F, Z directions. (In the case of an airplane Y is in tin? direction, of the 
span of the wings, Fis in the direction of flight, Z is vertical.) 


<7 = 

b 

t 

h 

f 
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W 




= dynamical pressure. 


=span of a wing (“ Broite”). 

— chord of a wing (“Tiefe”). 

=gap of a biplane (“ Hohe”). 

=area of surface (— h . t) (“Fliichc”). 
= lift (“ Auftriel) ”). 

=drag (“ Widerstand”). 


c w 


A 

Fq 

w 

Fq 


= lift. coefficient (—2 K y * 
— drag coefficient (=2 K x 


‘absolute”), 
“absolute ”). 


a = angle of attack. 


r = circulation. 

—velocity potential. 
=stream function. 
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